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^] . ABSTRACT 

' We present new optical and near infrared (NIR) photometry and spectroscopy of the type IIP 

supernova (SN), SN 2004et. In combination with already published data, this provides one of 
the most complete studies of optical and NIR data for any type IIP SN from just after explo- 
sion to +500 days. The contribution of the NIR flux to the bolometric light curve is estimated 
to increase from 15 per cent at explosion to around 50 per cent at the end of the plateau and 
then declines to 40 per cent at 300 days. SN 2004et is one of the most luminous IIP SNe which 
has been well studied and characterised, and with a luminosity of log L = 42.3 erg s and a 
56 Ni mass of 0.06 ± 0.04 M Q , it is 2 times brighter than SN 1999em. We provide parametrised 
bolometric corrections as a function of time since explosion for SN 2004et and three other 
IIP SNe that have extensive optical and NIR data. These can be used as templates for future 
events in optical and NIR surveys without full wavelength coverage. We compare the physi- 
cal parameters of SN 2004et with those of other well studied IIP SNe and find that the kinetic 
energies span a range of 10 50 -10 51 ergs. We compare the ejected masses calculated from hy- 
drodynamic models with the progenitor masses and limits derived from prediscovery images. 
Some of the ejected mass estimates are significantly higher than the progenitor mass esti- 
mates, with SN 2004et showing perhaps the most serious mass discrepancy. With the current 
models, it appears difficult to reconcile 100 day plateau lengths and high expansion velocities 
with the low ejected masses of 5-6 M implied from 7-8 M progenitors. The nebular phase 
is studied using very late time HST photometry, along with optical and NIR spectroscopy. The 
light curve shows a clear flattening at 600 days in the optical and the NIR, which is likely due 
to the ejecta impacting on circumstellar material. We further show that the [Oi] 6300, 6364 
A line strengths in the nebular spectra of four type IIP SNe imply ejected oxygen masses of 
0.5-1.5 M . 
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1 INTRODUCTION 

Supernovae (SNe) are classified based on the elements present in 
their spectra, with type II SNe being disting uished from Type I SNe 
by the presence of hydrogen (for review see lFilippenkdl 19970 . The 
type of explosion is directly determined by the evolutionary status 
of the progenitor star, which is influenced by the initial mass , metal- 
licity , rotation rate and presence of a binary companion (Smartt 
2009). The progeni tor stars of some type II SNe have b een shown to 
be red supergiants (Smartt et al. 2004; Li et al. 2005) and in some 
rare c ases blue supergiants l Walborn et al. 1987j: IPastorello et al.l 
2005). Stellar evolutionary models predict that most single stars 
with masses in the range ~ 8-30 M Q sh ould end their lives as 
red supergiants and produce type II SNe ([E ldridge & Tout 2004 ; 
lHegeretalJ|200i l Hirsc hi et al. | |2004 iLimongi & Chieffil 12003). 
This is supported by the observational constraints on core-collapse 
supernova progenitors, although the la ck of the detection of high 
mass progenitors is concerning ( Smartt 2009). 

Type II SNe can be further subclassified based on the shape of 
their light curves. Those that display a linear decrease from peak 
magnitude are called type IIL SNe, while those that display an ex- 
tended plateau of nearly constant luminosity are termed type IIP 
SNe. The plateau phase generally lasts ~ 80-120 days before en- 
tering the exponential decay phase and is caused by the diffusion 
of thermal energy deposited by the shock wave and by the release 
of internal energy when hydrogen starts to recombine. During re- 
combination, constant luminosity is achieved by the balance of the 
increase in radius caused by expansion and the inward movement 
(in Lagrangian coordinates) of the recombination front. It is obser- 
vationally established that the majority of IIP SN progenitors are 
red supergiants that initially have masses above 8 M Q and up to 
at least 17 M^. I Smartt et alj ( 120091) summarised the observed pro- 
genitor properties of the nearest IIP SNe and found a lack of high 
mass red supergiants in this sample. If these estimates of progenitor 
masses are accurate, it means long plateau phases are produced by 
ejected envelope masses as low as 6-7 M Q , which are in some cases 
much lower values than those c alculated from light curve modelling 
faamuy||200l : lNadezhinll2003T) . 

Type IIP SNe have been proposed as cosmological standard 
ca ndles and their u s e as di st ance indicator s has b e en demonstrated 
in IHamuv & Pintol d2002l) ; iNugent et all d2006l) ; IPoznanski et~ai1 
(2009). Therefore a detailed understanding of their pre-explosion 
parameters and explosion physics is of the utmost importance. 
Despite being the most common type of SN observed, there is 
a surprisingly small sample of IIP SNe to date with well mon- 
itored optical light curves and spectra from just after explosion 
through to the radioactive tail phase. In the past twenty years, 
there have been only four IIP SNe for which extended photo- 
metric and spectral monitoring observations have be en published; 
SN 1990E (fSchmidt et al] | 19931; iBen etti et alJI 19941), SN 1999em 
( lElmhamdi et alj|2003l; [Leonard et alj|2002al ; IHamuv et alj|200lh. 



SN 1999gi ^Leonard et al.ll2002bl) and SN 2005cs jPastorello etahl 
20061; iTsvetkov et al.l|2006; Bro wn et al.ll2007l ; lDessart et alj|2008l ; 



Pastorello et al. 2009). Coverage at near infrared wavelengths is 



even less common, with only one type IIP SN having been ob- 
served extensively both photometrically and spectr oscopically at 
these wavelengths in the last ten years, SN 1999em dHamuv et all 
200 ll ; lElmhamdi et alJl2003l ; lKrisciunas et alj2009l) . 

SN 2004et was discovered in the nearby starburst 
galaxy, NGC 6946 by S. Mo retti on 2004 September 27 
dZ witter. Munari & Morettil |2004|). At a distance of on ly 5.9 ± 
0.4 Mpc l Karachentsev. Sharina & Huchtmeieil |2000|) . it was 



an ideal candidate for intensive follow-up observations. A high 
resolution spectrum of SN 2004et was obtained with the Mt. 
Ekar 1.82-m telescope on 2004 September 28 that confirmed it 
was a type II event with sp ectra showing prominent H B ai mer 
lines with P-Cygni profiles dZwitter. Munari & Morettil l2004h . A 
total extinction of E(B-V) = 0.4 1 ± 0.007 was also estimated by 
IZwitter. Munari & Morettil |2004) from the equivalent width of the 
NalD lines in the high resolution spectrum. SN 2004et (R.A.: 20 h 
35 m 25 s . 33, Deck: +60° 07' 17".7 (12000)) was located in one of 
the spiral arms of the galaxy, which are known to be regions of 
high star formation. NGC 6946 is presently the galaxy with the 
highest number of SNe discovered to date, 9 SNe since the first in 
SN 1917A with the most recent be ing SN 2008S, which has been 
suggested bv lBotticella et alj d2009h as a pr obable electron-ca pture 
even t, although i ts natu re has been debated dSmifh et alj2009l) . 

ISahu et alj d2006l) reported optical photometric and spectro- 
scopic observations of SN 2004et from 8-541 days post explosion. 
They noted that SN 2004et was at the brighter end of SNe IIP lu- 
minosities and estimated the mass of 56 Ni synthesised during the 
explosion to be 0.06 + 0.02 M Q . They also suggested that the steep- 
ening of the decline rates of the optical luminosity one year after 
explosion along with a blueshift in the e mission lines at a s imilar 
epoch was an indication of dust formation. iMisra et alj d2007h pre- 
sented optical photometry of SN 2004et from ~ 1 4^470 days post 
explosion and in agreement with lSahu et al. I d2006l) . estimated an 
ejected 56 Ni mass of 0.06 ± 0.03 M Q . 

SN 20 04et was detected at r adio frequencies just 14 days post 
explosion dStockdale etai1l2004h and this early detection suggests 
the pres ence of a ppreciable circumstellar material (CSM) around 
the SN dSahu et alj|2006f) . The SN was extensively monitored by 
IStockdale etall d2004l) u sing the VLA at 22.4 and 8.4 GHz and 
by iBeswick et alj d2004l) using MERLIN at 4.9 GHz. SN 2004et 
was a lso observed on 2 005 lanuary 02 using the GMRT at 1.4 
GHz (Mis ra et al ]|2007l) and an 8.4 GHz VLBI observation was 
obtained by iMartf-Vidal et alj ( 120071) on 2005 Feburary 20 show- 
ing a clear asymmetry in the emission structure that can be ex- 
plained if the SN ejecta expanded in a clumpy CSM. The radio 
lumi nosity of SN 2004et was found to be am ong the highest for IIP 
SNe dChevalier. Fransson & Nvm ark 2006) with an estimated pre- 
supernova mass-loss rate of 1.5-3 x 10~ 6 M Q yr~' (using a wind 
velocity of 10 km s~'), which suggests a progenitor mass of ~ 20 
M . X-ray data of SN 2004et were also obtained using th e Chandra 
X-ray Observatory at 30 , 45 and 72 days post explosion dRho et all 
120071 ; IMisra et al]|2007l) . The unabsorbed X-ray luminosity sug- 
gests a wind density parameter value of 2.5, which is a factor of 
2 larger than the wind dens ities calculated for SN 1999em and SN 
2004dj dChugai et alj|200 7). This gives a pre-supernova mass-loss 
rate of 2-2.5 x 10~ 6 M Q yr~', which is consistent with estimate 
from the radio observa tions. 

Kota k et al. I d2009h presented mid infrared (MIR) observations 
of SN 2004et obtained with the Spitzer Space Telescope ranging 
from 64-1240 days post explosion, along with three very late- time 
optical spectra. They reported spectroscopic evidence for silicate 
dust formation in the ejecta of SN 2004et with a total mass for the 
dust of mass of a few times 10~ 4 M Q , which would not make a ma- 
jor contribution to the total mass of cosmic dust. The most promi- 
nent spectral emission lines in the very late-time optical spectra are 
observed to display boxy profiles, which was suggested to be a sig- 
natur e of ejec ta-CSM interaction. 

ILi et alj d2005l) reported that the progenitor of SN 2004et had 
been identified in a pre-expl osion optical image o f NGC 6946 as a 
yellow supergiant. However. ICrockett et alj (2009) have shown that 
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the source indicated as the progenitor bv lLietal.l ( l2005h is still vis- 
ible in images taken after SN 200 4et has faded and so c annot be the 
progenitor star. At the same time, ICrockett et all [2009) also identi- 
fied an alternative progenitor star in a pre-explosion i band image 
taken at the Isaac Newton Telescope with the Wide Field Camera. 
This star was not detected in pre-explosion V or R band images 
suggesting that the progenitor of SN 2004et is either intrinsically 
red or surrounded by dust. The luminosity of the progenitor was 
estimated using the / band photometry and R-I colour limits. Us - 
ing the STARS stellar evolutionary model feldridge & Toutll2004l) . 
the progenitor star was found consistent with a lat e K to late M- 
type s upergiant with an initial mass of 8+j M Q (see lCrockett et all 
(2009) for more details). 

In this paper, we present our own complete data set for SN 
2004et making one of the most comprehensive sets of optical and 
NIR photometry and spectroscopy for a type IIP SNe, with data 
coverage from soon after explosion through to the radioactive tail 
phase. The information obtained from analysis of the light curve 
and spectra is compared to the progenitor information obtained 
from the pre-explosion images to explore if the resulting progen- 
itor properties are consistent. Extensive NIR spectroscopy is rare 
for type IIP SNe and so we investigate what further results can be 
obtained from analysis of NIR spectral features. SN 2004et is also 
included in a wider sample of nearby type IIP SNe to determine 
how its properties fit in the overall picture for the core-collapse of 
massive, hydrogen-rich stars. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Optical photometry 

Optical photometric observations of SN 2004et were carried out us- 
ing a range of telescopes, beginning soon after its discovery (2004 
September 27) until over three years post explosion (2008 January 
18). Observations at very early times were taken with a 0.41-m tele- 
scope equipped with an Apogee AP47p CCD in BVRI filters and 
the 1.82-m Copernico telescope equipped with AFOSC, located at 
Mt. Ekar-Asiago (Italy). Three epochs of data were taken with the 
3.58-m Telescopio Nazionale Galileo (TNG) using DOLORES, lo- 
cated at Roque de Los Muchachos Observatory, La Palma, Canary 
Islands, along with data from the 0.72-m TNT at the Collurania Ob- 
servatory (Teramo, Italy). Photometric observations were also ob- 
tained with the 0.7-m AZT2 telescope at the Sternberg Astronom- 
ical Institute, Moscow, the 0.6-m Z600 telescope at the Crimean 
Laboratory of SAI, Nauchny, Crimea, the 0.38-m KGB telescope 
at the Crimean Astrophysical Observatory, Nauchny, Crimea and 
the 1-m SAO Z1000 telescope at the Special Astrophysical Obser- 
vatory of the Russian Academy of Sciences, Zelenchuk. Two late 
time observations were also taken with the 2.56-m Nordic Opti- 
cal Telescope (NOT) at Roque de Los Muchachos Observatory, La 
Palma, Canary Islands using ALFOSC. 

All the images were trimmed, bias subtracted and flat-field 
corrected using the standard iraiQ tasks. A sequence of stars in the 
field of view of SN 2004et were calibrated with re spect to the se- 
quence star magnitudes given in lMisra et ai]( l2007l) . The sequence 



iraf is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under the cooperative agreement with the National Science 
Foundation. 




Figure 1. The field of SN 2004et with the supernova and comparison stars 
marked. This image was taken with the 1.82m telescope at Mt. Ekar, Italy 
~ 3 weeks post explosion in the V band. 



star magnitudes of iMisra et al.l d2007t) were also confirmed inde- 
pendently using standard star observations obtained on some of the 
nights. The errors on the reported magnitudes for SN 2004et were 
determined by adding in quadrature the errors on the instrumental 
magnitudes and the errors due to the magnitude calibration. Figure 
Q]shows a finding chart for SN 2004et, with the secondary standard 
stars marked. 

Additional late time Hubble Space Telescope (HST) WFPC2 
observations were obtained as part of the proprietary program GO- 
11229 (PI: Meixner) in July 2007 and January 2008, which have 
since been made public. These images were downloaded from the 
Space Telescope Science Institute (STScI) archive using the OTFR 
pipeline. Observations were obtained in two filters F606W and 
F814W that are equivalent to the Johnson V and / bands respec- 
tively. PSF-fitting p hotometry was performed on the data using the 
HSTphot package (Dolphin 2000). HSTphot corrects for chip-to- 
chip variations and aperture corrections. No colour corrections have 
been applied to convert to V and / filters because the colour trans- 
formations to the Johnson-Cousins magnitude system are not well 
constrained for late time SN spectra. The magnitude values quoted 
in Table[T]are for the HST filters in the Vegamag system. One epoch 
of very late time B VRI data was also obtained with the William Her- 
schel Telescope (WHT) at Roque de Los Muchachos Observatory, 
La Palma, Canary Islands using the Auxiliary Port Imager (AUX). 

2.2 Near infrared photometry 

SN 2004et was observed in the NIR JHK bands with SWIRCAM 
mounted at the focal plane of the 1.08-m AZT-24 telescope, at 
Campo Imperatore Observatory (Italy) on 27 epochs between 2004 
September 30 and 2006 September 21. The AZT-24 telescope is op- 
erated jointly by the Pulkovo observatory (St. Petersburg, Russia) 
and INAF Observatorio Astronomico di Roma/Collurania (Campo 
Imperatore, Italy). SWIRCAM is equipped with a Rockwell PIC- 
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Table 1. Log of optical photometric observations of SN 2()()4et. 



Date 


JD (2450000+) 


Phase* (days) 


U 


B 


V 


R 


/ 


Source 


o r\ /nn maa a 

29/09/2004 


53278.0 


7.5 




I "i c\ A , A Al 

12.94±0.01 


1 o, s i , A A1 

12.67±0.01 


1 O O C .A AO 

12.35+0.02 


1 o 1 o . n no 

12.13+0.02 


TMT A 1 

IN 1 0.72m 


01/10/2004 


53279.9 


9.4 




12.90±0.04 


12.60±0.02 


1 O O O . A AO 

12.28+0.02 


1 O A 1 . A AO 

12.01+0.02 


0.41m+AP47p 


i\ a 1 a a/oaa,i 

04/10/2004 


53283.0 


12.5 




a ^ r\o , n f\A 

12.98±0.01 


1 O £1 C i A AO 

12.65±0.02 


1 O T7 i A A 1 

12.27+0.01 


1 O A 1 i A A 1 

12.01+0.01 


T'TVTO - ' n TO 

INT 0.72m 


rx^i / 1 a /o aa a 

07/10/2004 


C O O O /*"" f\ 

53286.0 


15.5 




13.02±0.01 


1 O /I A i A A 1 

12.64±0.01 


1 o o /T i A no 

12.26+0.02 


i 1 no i n no 

1 1.98+0.02 


TT^TT' A 7A 

INT 0.72m 


1 A / 1 A /O A A A 

10/10/2004 


c o o o o n 

53288.9 


18.4 




13.00±0.02 


1 O CO i A AO 

12.58±0.02 


1 O O O , A A 1 

12.23+0.01 


i 1 no i n no 

1 1.98+0.02 


f\ A A A A 1 

0.41m+AP47p 


1 r /1 a /o a a a 

15/10/2004 


53294.2 


23.7 


13.04±0.04 


13.16+0.01 


1 O f A . A A 1 

12.61+0.01 


12.16+0.01 


11 A A . A A 1 

ll. 94+0.0 1 


AZT2 0.7m 


1 "7 1 A A /O A A A 

17/10/2004 


53295.9 


25.4 




13.27±0.02 


1 O y^O . A A 1 

12.63±0.01 


1 O O O . A A 1 

12.22+0.01 


1 1 AO . A A 1 

H.92+0.01 


0.4lm+AP47p 


19/10/2004 


53297. 1 


26.5 








1 O O A . A A 1 

12.24+0.01 


11 A A . A A 1 

H.94+0.01 


Ekar 1.82m+AFOSC 


i n m n /o a a a 

19/10/2004 


^ t AAn n 

53298.0 


27.5 




13.41±0.02 


-1 O y — J . A A -f 

12.67±0.01 


1 O O O . A A 1 

12.22+0.01 


11 o y . A A 1 

H.86+0.01 


TNT 0.72m 


OO IA A /O AA A 

28/10/2004 


53307.4 


36.9 


i a •"> "7 . n r\ A 

14.37±0.04 


13.79±0.01 


1 O 1 A . A. A 1 

12.79±0.01 


i o o /T . n n 1 

12.26+0.01 


11 (~\ C , A A 1 

1 1.95+0.01 


AZ12 0.7m 


Til/1 A /O A A A 

30/10/2004 


53309.4 


38.9 


14.51±0.03 


13.86+0.01 


12.81+0.01 


12.29+0.01 


11.95+0.01 


Z600 0.6m 


02/1 1/2004 


53312.4 


41.9 


14.70±0.09 


13.93±0.01 


12.83+0.01 


12.27+0.01 


11.95+0.01 


AZT2 0.7m 


05/1 1/2004 


53315.4 


54.9 


14.85±0.10 


14.04±0.01 


12.87±0.01 


12.30+0.01 


11.97+0.01 


Z600 0.6m 


06/1 1/2004 


53316.3 


45.8 


15.08±0.03 


14.10±0.01 


12.95±0.01 


12.40+0.01 


12.01+0.01 


Z600 0.6m 


A"7 fl 1 /a/»n A 

07/11/2004 


COO 1 T O 

53317.2 


46.7 


15.07 ±0.01 


A A f\£> , f\ f\A 

14.08±0.01 


1 O O O , A. A 1 

12.88±0.01 


1 o oo . n n 1 

12.32+0.01 


11 A C , A A 1 

11.95+0.01 


'■"z/rnn n f 

Z600 0.6m 


AO M i /o a a A 

08/1 1/2004 


c o o i ~j a 

53317.9 


47.4 




1 A 1 T , A At 

14.12±0.01 


1 O O A . A A 1 

12.89±0.01 


1 O O /~ . A A 1 

12.36+0.01 


11 AO . A A 1 

11.93+0.01 


f\ A A A A -7 

0.41m+AP47p 


a o I a i /o a a a 

08/1 1/2004 


53318.5 


48.0 


15.09 ±0.17 


I4.08±0.0l 


1 o a A . A A 1 

12.89±0.01 


1 O O O . A A 1 

12.38+0.01 


1 1 . A A 1 

H.97+0.01 


Z600 0.6m 


a a 1 1 i /o a a 

09/1 1/2004 


53319.4 


48.9 


AC A C . /\ /\1 

15.15±0.07 


1 A 1 ^> . A Al 

I4.13±0.0l 


1 O A 1 . A A 1 

12.91+0.01 


1 O O A . A A 1 

12.34+0.01 


11 A S~ . A A 1 

ll. 96+0.0 1 


Z600 0.6m 


1AM A /OA A A 

10/11/2004 


j— -) o O A a 

53320.0 


49.5 






1 O r\C . A AO 

12. 95+0.03 


12.37+0.03 




TNT 0.72m 


10/11/2004 


53320.4 


49.9 


15. 25 ±0.09 


14.15+0.01 


12.90±0.01 


12.34+0.01 


11 f\C . A A 1 

H.95+0.01 


Z600 0.6m 


14/1 1/2004 


53323.9 


52.6 


15.38±0.02 


I4.25±0.0l 


12.91±0.02 


12.36+0.01 


H.91+0.01 


TNG+LRS 


16/11/2004 


53325.8 


55.3 




I4.30±0.0l 


12.99 ±0.01 


12.39+0.01 


H.92+0.01 


Ekar 1.82m+AFOSC 


16/11/2004 


53326.0 


55.5 




14.22±0.02 


12.94±0.03 


12.36+0.02 


ll. 92+0.03 


TNT 0.72m 


1 O 1 A A /OA A A 

18/11/2004 


j— -) o O o a 

53328.0 


57.5 




A A TA . A A^ 

14.30±0.02 


1 O A A . A AO 

12.99±0.02 


1 O O O . A A 1 

12.38+0.01 


11 A O . A A 1 

11.93+0.01 


TNT 0.72m 


A C\ 1 A A /O A A A 

19/11/2004 


Co -> o A / \ 

53329.0 


58.5 




14.36±0.02 


-f O A -1 . A A 1 

13.01+0.01 


1 O A A . A A 1 

12.40+0.01 


11 A 1 . A A 1 

H.91+0.01 


T^l 1 OO A "T^/~\0^1 

Ekar 1.82m+AFOSC 


a r\ I a a mt~\f\ a 

19/11/2004 


-> -) O A o 

53329.2 


58.7 


15.65±0.02 


14.32±0.01 


1 O A ~T . A A 1 

12.97+0.01 


1 O O . A A 1 

12.36+0.01 


11 A y" . A A 1 

H.96+0.01 


Z600 0.6m 


OA 1 A A /OA A A 

20/1 1/2004 


j— -> o o a o 

53329.8 


59.3 




-\ A 1A , A Al 

14.39±0.01 


1 O AO . A A 1 

12.98+0.01 


1 O O O . A A 1 

12.38+0.01 


11 A O . A A 1 

H.93+0.01 


n i 1 a "n a -7 

0.4lm+AP47p 


n 1 / 1 i /o a a a 

21/11/2004 


j— -) o O 1 O 

53331.2 


60.7 


15.60±0.07 


14.37±0.02 


1 O AO . A A 1 

12.98+0.01 


1 O 1 . A A 1 

12.37+0.01 


1 1 A ~7 . A A 1 

11.97+0.01 


Z600 0.6m 


22/1 1/2004 


53332.0 


61.5 




14.40±0.02 


13.02+0.03 






TNT 0.72m 


OO 1 A A /OA A A 

22/1 1/2004 


53332.2 


61.7 




14.45±0.04 


12.95+0.03 


12.31+0.02 


1 1 AO . A AO 

ll. 92+0.03 


KGB 0.38m 


23/1 1/2004 


53333.0 


62.5 




14.38±0.02 


13.02+0.03 


12.42+0.02 


11.91+0.02 


TNT 0.72m 


25/1 1/2004 


53335.3 


64.8 


A C f\ CI .1X1^ 

I5.98±0.l2 


14.42±0.02 


13.01+0.01 


12.38+0.01 


11.97+0.01 


AZT2 0.7m 


1 A / 1 O /OA A A 

10/12/2004 


53349.8 


79.3 




14.73±0.01 


13.14+0.01 


1 O AC . A A 1 

12.45+0.01 


11.93+0.01 


T^l 1 OO A "T^/~\0^1 

Ekar 1.82m+AFOSC 


1 O M O /OAA A 

13/ 12/2004 


f n r a *7 

53352.7 


82.2 




1 A n . c\ t~\ A 

14.77±0.0l 


1 O 1 A . A A 1 

13.14+0.01 


1 O ,1-7 . A A 1 

12.47+0.01 




T^l 1 OO A "T^/~\0^1 

Ekar 1.82m+AFOSC 


1 O / 1 O /OA A A 

13/12/2004 


f T T r o A 

53353.0 


82.5 




14.68±0.01 


1 O 1 y . A A 1 

13.16+0.01 


1 O A A . A A 1 

12.49+0.01 


1 1 A ~7 . A AO 

11.97 ±0.02 


TNT 0.72m 


14/12/2004 


53354.0 


83.5 




14.68±0.01 


13.15+0.01 


12.52+0.02 


11.99+0.01 


TNT 0.72m 


1 £■ / 1 O /OA A A 

15/12/2004 


53354.8 


84. 3 




A A OO * A Ad 

14.88±0.04 


1 O IT . A AO 

13.17+0.03 


1 O A A . A AO 

12.49+0.03 


11 AO . A AO 

11.93 ±0.03 


T^l 1 OO A "T^/~\0^1 

Ekar 1.82m+AFOSC 


23/12/2004 


53363.0 


92.5 




14.86±0.02 


13.24+0.02 


12.57+0.02 


11.99±0.02 


TNT 0.72m 


28/12/2004 


coo/ t > a 

53368.0 


97.5 




1 A AO . A AO 

14.93±0.03 


13.28+0.01 


1 O f A . A A 1 

12.61+0.01 


12.10±0.01 


TNT 0.72m 


29/12/2004 


Coo x ( \ -> 

53369.3 


98.8 






13.32+0.01 


12.56+0.01 


12.12+0.01 


AZT2 0.7m 


A -> tr\ \ ir\ n A c 

03/01/2005 


53373.8 


103.3 




AC A A i C\ f\ A 

15.14±0.04 


1 O A "1 i A A 1 

13.41+0.01 


1 O /T A i A (\ A 

12.69+0.01 


1 O 1 /T i A A 1 

12. 16±0.01 


7 ' 1 1 OO i A T^/^v O 

Ekar 1.82m+A10SC 


I\C tf\ 1 /OAAC 

05/01/2005 


53376.0 


105.5 




AC 1 O ■ f\ f\ A 

15.18±0.04 


i O A f . A AO 

13.46+0.03 


1 o oo . n n A 

12.83+0.04 


1 O 1 n , A A£T 

12.10+0.05 


TNT 0.72m 


t\i /a i /oaac 

07/01/2005 


c o o "7 o a 

53378.0 


107.5 




1 C O O , A AO 

15.23±0.02 


1 O C A , A AO 

13. 54+0.02 


1 O OA . A AO 

12.80+0.02 


1 O OT . A AO 

12.27±0.03 


TNT 0.72m 


1 1 /A A /O A A C 

11/01/2005 


C O o o o o 

53382.2 


111.7 




AC A 1 . A A 1 

15.47±0.01 


1 O — T A . A A 1 

13.70+0.01 


1 O O O . A A 1 

12.82+0.01 




AZT2 0.7m 


14/01/2005 


COO O ,1 O 

53384.8 


1 14.3 




AC SO . A A 1 

15.68±0.01 


i O O i , A A 1 

13.81+0.01 


i o no . n n 1 

13.02+0.01 


i o o "7 . n n i 

12.37±0.01 


Ekar l.82m+AEOSC 


1 "7 /A 1 ^MAr 

17/01/2005 


c o o o o a 

53388.0 


117.5 




1 C "1 A .A AO 

15.74±0.03 


1 O A C .A AO 

13. 95+0.03 


1 O 1 A . A AO 

13.19+0.03 


1 O C /~ . A AO 

12.56±0.03 


TNT 0.72m 


1 n /A 1 /1AAC 

18/01/2005 


coo on 1 

53389.3 


118.8 




1 c in . a ac 

15.79±0.05 


1 A AA . A AO 

14.09+0.03 


1 o ii , A no 

13.1 1+0.02 


12.64±0.02 


A TTO A. ^7 

AZ12 0.7m 


f\ lf\ A /OA AC 

20/01/2005 


53391.5 


120.5 




1 f A f . A A A 

16.06±0.04 


14.21+0.04 


13.42+0.04 


12.77±0.04 


TNT 0.72m 


A A /AO ^AAr 

04/02/2005 


53406.2 


135.7 




17.34±0.05 


15.64+0.03 


14.47+0.03 


13.86±0.02 


AZT2 0.7m 


1 A /AO /OAAC 

10/02/2005 


53412.0 


141.5 






15.66+0.06 






TNI 0.72m 


10/02/ZOOj 


53412.2 


141.7 




1 /.3o±o.oy 


1 J. /j+O.OD 


1 A CA_lO AO 

t4. JO+O.Oz 


1 1 m_i_o m 
t3.y /±0.02 


A 7T1 O Tm 

AZtz U. /m 


03/03/2005 


53432.5 


162.0 




17.63±0.07 


15.93+0.04 


14.67+0.03 


14.14±0.02 


AZT2 0.7m 


09/03/2005 


53439.2 


168.7 




17.74±0.02 


15.97+0.01 


14.89+0.01 


14.06±0.01 


Ekar 1.82m+AFOSC 


17/03/2005 


53446.6 


176.1 




17.55±0.06 


16.12+0.04 


14.81+0.03 


14.26±0.03 


AZT2 0.7m 


18/03/2005 


53447.0 


176.5 




17.75±0.01 


16.04+0.01 


14.96+0.01 




Ekar 1.82m+AFOSC 


27/03/2005 


53456.5 


186.0 




17.77±0.05 


16.14+0.03 


14.92+0.03 


14.39±0.02 


AZT2 0.7m 


29/03/2005 


53458.5 


188.0 




17.77±0.08 


16.21+0.05 


14.92+0.03 


14.41±0.02 


AZT2 0.7m 


08/04/2005 


53468.5 


198.0 




17.78±0.17 


16.24+0.12 


14.96+0.05 


14.47±0.06 


AZT2 0.7m 


13/04/2005 


53473.5 


203.0 




17.87±0.05 


16.35+0.03 


15.06+0.03 


14.53±0.03 


AZT2 0.7m 


18/05/2005 


53509.5 


239.0 




17.94±0.07 


16.69+0.04 


15.46+0.03 


14.99±0.03 


AZT2 0.7m 


15/06/2005 


53537.0 


266.5 




18.35±0.02 


16.97+0.02 


15.76+0.01 


14.99±0.03 


Ekar 1.82m+AFOSC 


02/07/2005 


53554.0 


283.5 




18.44±0.03 


17.13+0.02 


15.88+0.01 


15.24±0.02 


Ekar 1.82m+AFOSC 


13/07/2005 


53565.4 


294.9 




18.48±0.07 


17.39+0.04 


16.11+0.03 


15.74±0.03 


AZT2 0.7m 



Optical and near infrared coverage of SN 2004et 5 



Table 1 - continued 



Date 


JD (2450000+) 


Phase* (days) 


V 


B 


V 


R 


/ 


Source 


14/07/2005 


53566.5 


296.0 




18.52±0.03 


17.20±0.03 


16.14±0.03 


15.83±0.03 


SAO Z1000 lm 


28/07/2005 


53580.4 


309.9 




18.62±0.08 


17.33±0.03 


16.27±0.03 


15.99+0.12 


AZT2 0.7m 


04/08/2005 


53587.4 


316.9 




18.70±0.06 


17.53±0.04 


16.34±0.03 


16.08±0.03 


AZT2 0.7m 


10/08/2005 


53593.4 


322.9 




18.75±0.06 


17.47±0.05 


16.40±0.03 


16.19±0.03 


AZT2 0.7m 


16/08/2005 


53599.4 


328.9 




18.86±0.07 


17.59±0.10 


16.49±0.03 


16.23±0.04 


AZT2 0.7m 


24/08/2005 


53607.4 


336.9 




18.84±0.13 


17.76±0.08 


16.65±0.03 


16.42±0.03 


AZT2 0.7m 


28/08/2005 


53611.2 


340.7 


20.17±0.03 


18.91±0.02 


17.78±0.02 


16.71+0.01 


16.12+0.01 


TNG+LRS 


30/08/2005 


53613.4 


342.9 




18.84±0.05 


17.77±0.03 


16.68±0.03 


16.53±0.05 


Z600 0.6m 


31/08/2005 


53614.5 


344.0 




18.67±0.07 


17.87±0.05 




16.43±0.09 


Z600 0.6m 


07/09/2005 


53621.3 


350.8 




18.8O±0.O8 


17.95±0.12 


16.76±0.04 


16.66+0.06 


AZT2 0.7m 


14/09/2005 


53628.1 


357.6 




19.08±0.03 


17.97±0.02 


17.03±0.02 


16.39+0.02 


NOT+ALFOSC 


14/09/2005 


53628.4 


357.9 




18.95±0.06 


17.84±0.06 


16.81±0.04 


16.61+0.05 


AZT2 0.7m 


03/10/2005 


53647.4 


376.9 






18.32=t0.18 


17.15±0.05 




AZT2 0.7m 


10/10/2005 


53654.0 


383.5 




19.25±0.02 


18.24±0.01 


17.39±0.02 


16.71+0.09 


Ekar 1.82m+AFOSC 


25/10/2005 


53669.0 


398.5 




19.40±0.02 


18.42±0.01 


17.48±0.01 


16.76+0.02 


Ekar 1.82m+AFOSC 


01/11/2005 


53676.4 


405.9 




19.44±0.06 


18.54±0.04 


17.54±0.03 


17.32+0.13 


Z600 0.6m 


03/11/2005 


53677.8 


407.3 




19.47±0.03 


18.56±0.01 


17.56±0.01 


16.86+0.04 


Ekar 1.82m+AFOSC 


07/11/2005 


53682.3 


411.8 




19.45±0.05 


18.60±0.04 


17.63±0.03 


17.62+0.14 


Z600 0.6m 


09/11/2005 


53684.3 


413.8 




19.27±0.09 


18.62±0.10 


17.71±0.04 




Z600 0.6m 


21/05/2006 


53876.2 


605.7 






21.43±0.04 


20.81±0.02 




NOT+ALFOSC 


30/09/2006 


54008.0 


737.5 






22.45±0.08 






NOT+ALFOSC 


08/07/2007 


54289.5 


1019.0 






23.01+0.01" 




22.85+0.02" 


HST+WFPC2 


12/08/2007 


54324.5 


1054.0 




23.62±0.06 


23.06±0.04 


22.27±0.05 


21.99+0.06 


WHT+AUX 


18/01/2008 


54484.8 


1214.3 






23.24+0.01° 




22.13+0.02° 


HST+WFPC2 



* relative to the epoch of date of explosion (JD = 2453270.5) 
°HST (F606W = V band, F814W =* / band) 



NIC array with 256 x 256 pixels with a pixel size of 1.04". One 
epoch of JHK data was also obtained on 2005 July 24 at the 3.58- 
m Telescopio Nazionale Galileo (TNG), using NICS with a 1024 x 
1024 pixel Rockwell array with a pixel size of 0.25". At each epoch, 
multiple dithered images of the SN were acquired along with mul- 
tiple dithered sky images. The sky images were median combined 
to eliminate stars. The resulting sky frames were subtracted from 
the SN frames for each band. The sky-subtracted SN images were 
then aligned and co-added. The data were processed using standard 
tasks in iraf. 

The photometric calibration was carried out using the magni- 
tudes of nearby stars taken from the Two Micron All Sky Survey 
(2MASS) photometric catalogue. Point spread function (PSF) fit- 
ting magnitude measurements were performed using iraf's daophot 
package. Differential photometry of the SN magnitude was then 
carried out by comparison with the calibration of the 2MASS stan- 
dards. The 2MASS filters differ from those of SWIRCAM and 
NICS, particularly in the K s bands but the colour terms derived 
from the reference stars were small, which implies only a negli- 
gible difference in photometric systems. The resulting SN magni- 
tudes and associated errors are given in Table|2] 

Late time NIR data had been obtained with HST using 
NICMOS in July 2007 and January 2008 as part of the propri- 
etary program GO-11229 (PI: Meixner). Multiple exposures were 
taken in each of the filters, F110W, F160W and F205W, that had 
been offset by sub-pix el amounts, which could then be 'drizzled' 
( Fruchter & Hook 2002) to improve the spatial sampling of the 
PSF. PSF-fitting photometry was performed on the NIR 'drizzled' 
data using iraf's daophot package. The filters used, F110W, F160W 
and F205W are roughly equivalent to the /, H and K filters respec- 
tively. However no conversion equations have been applied to con- 
vert to J, H and K filters since the colour transformations are not 



well constrained for late time SN spectra (in fact no NIR spectra 
at this epoch exist for accurate S-correction). The values quoted in 
Table|2]are for the HST filters in the Vegamag system. 



3 PHOTOMETRIC EVOLUTION 
3.1 Light curve 

Figure|2]shows the optical and NIR light curve (UBVRIJHK) of SN 
2004et from just after explosion to +414 days post explosion. Very 
early tim e R band magnitudes w ere reported by Klotz and collab- 
orators in Yama oka et alj ( 120041) . Nothing was visible to a limiting 
magnitude of 19.4 ± 1.2 at the SN position on 2004 September 
22.017 UT, when imaged by the robotic TAROT telescope. Pre- 
discovery detection of SN 2004et on September 22 .983 UT exists, 
with the SN having a magnitude of 15.17 ± 0.16 ( Yamao ka et al.l 
2004). The explosion epoch of SN 2004et has therefore been well- 
constrained and the explosion epoc h is taken as 2004 September 
22.0 (JD 2453270.5), the same as in lLi et alj fe005t) . 

The characteristic plateau of type IIP SNe is visible up to 1 10 
± 15 days, before showing a sharp decline onto the radioactive tail, 
whose luminosity is powered by the decay of 56 Co to 56 Fe. The ex- 
pected decline rate during the neb ular phase is 0.98 mag per 100 
days for complete y-ray trapping dPatat et all 1 1994b . The decline 
rates per 100 days for the BVRI bands were calculated using a least 
squares fit, during the early nebular phase (~ 136-300 days) and 
found to be y B = 0.66 ± 0.02, j v = 1.02 ± 0.01, j R = 0.92 ± 0.01 
and 7/ = 1.09 ± 0.01. These values are close to the expected de- 
cline rate and are consistent with the early nebular phase flux being 
dominated by radioactive decay. 

After 300 days, steeper decline rates were seen in the BVR 
bands and a comparison of the decline rates is shown in Figure [3] 
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Table 2. Near infrared photometric observations of SN 2004et. 



Date JD (2450000+) Phase* (days) J H K s 



30/09/2004 


53278.8 


8.3 


11.89±0.04 


11.58+0.05 


11.33+0.05 


01/10/2004 


53279.8 


9.3 


11.94±0.04 


11.55+0.06 


11.31+0.05 


02/10/2004 


53280.9 


10.4 


11.93+0.05 


11.54+0.06 




03/10/2004 


53281.9 


11.4 


11.85+0.04 


11.48+0.05 


11.32+0.05 


05/10/2004 


53283.9 


13.4 


11.79+0.04 


11.46+0.05 


11.27+0.05 


07/10/2004 


53285.8 


15.3 


11.76+0.04 


11.44+0.06 


11.19+0.06 


08/10/2004 


53286.8 


16.3 


11.74+0.04 


11.40+0.06 


11.16+0.06 


09/10/2004 


53287.7 


17.2 


11.81+0.04 


11.37+0.09 


11.03+0.08 


18/10/2004 


53296.9 


26.4 


11.64+0.04 


11.36+0.08 


10.92+0.08 


20/10/2004 


53298.8 


28.3 


11.59+0.05 


11.25+0.06 


10.97+0.07 


22/10/2004 


53300.8 


30.3 


11.60+0.04 


11.18+0.08 


10.96+0.06 


23/10/2004 


53301.8 


31.3 


11.54+0.05 


11.22+0.08 


11.03+0.06 


27/10/2004 


53305.8 


35.3 


11.55+0.03 


11.18+0.05 


10.96+0.05 


17/11/2004 


53326.8 


56.3 


11.42+0.05 


11.07+0.05 


10.93+0.05 


18/11/2004 


53327.9 


57.4 


11.44+0.04 


11.08+0.05 


10.88+0.05 


20/1 1/2004 


53329.7 


59.2 


11.44+0.04 


11.09+0.05 


10.88+0.05 


21/11/2004 


53330.7 


60.2 


11.43+0.04 


11.07+0.05 


10.84+0.05 


23/11/2004 


53332.8 


62.3 


11.48+0.04 


11.12+0.05 


10.89+0.05 


25/11/2004 


53334.8 


64.3 


11.45+0.04 


11.06+0.04 


10.87+0.05 


07/12/2004 


53346.7 


76.2 


11.47+0.05 


11.11+0.08 




13/12/2004 


53352.9 


82.4 


11.43+0.04 


11.08+0.05 


10.86+0.05 


15/12/2004 


53354.7 


84.2 


11.46+0.04 


11.08+0.05 




04/01/2005 


53374.7 


104.2 


11.64+0.04 


11.30+0.05 


11.08+0.05 


10/01/2005 


53380.8 


110.7 


11.69+0.05 


11.36+0.04 


11.16+0.05 


03/02/2005 


53404.7 


134.2 


13.21+0.04 


12.83+0.06 




10/02/2005 


53411.7 


141.2 


13.37+0.04 


13.06+0.07 




24/07/2005 


53576.6 


306.1 


15.36+0.14" 


15.69+0.12" 


15.50+0.16" 


21/09/2006 


54000.0 


729.5 


> 20.2 






08/07/2007 


54289.5 


1019.0 


22.08+0.02 6 


22.57+0.08 6 


> 22.3 6 


19/01/2008 


54485.6 


1215.1 


22.30+0.03'' 


22.70+0. 10 6 


21.18+0.06* 



*relative to the epoch of date of explosion (JD = 2453270.5) 
"TNG+NICS 

^HST (Fl 10W = / band, F160W =s H band, F205W =* K band) 



_' 1* 




• 


• U-4.5 • 1 + 3 
A B-2.5 * J+4.1 


: * Y 
v*. 


■ V • H + 4.9 


□ R+1.7 O K + 5.6 






- Zt & " 

% -- 


A 

• 




• 

• • 



200 300 400 

nee explosion (JD-2453270.5) 



The decline rates based on data from ~ 136-296 days only, are 
shown as a solid line for each of the bands and are extrapolated out 
to 800 days. The decline rates using data from between ~ 296-414 
days are y B = 0.85 ± 0.02, y v = 1.17 ± 0.02, y R = 1.33 ± 0.01 
and y; = 0.93 ± 0.02 and are shown as the dashed lines in Figure 
[3] Note that the late time / band photometry shows a large scatter 
due to the intrinsic differences in the instrumental / band filters and 
this could affect our estimates of the slope of the / band light curve. 
The points after 600 days shown in Figure [3] were not included in 
the determination of the slopes of the V and R bands. The BVR 
bands became steeper during this period while the slope of the / 
band became slightly less steep. This steepening could imply either 
the leakage of y rays because the SN had become transparent to y 
rays at this stage or dust formation in the SN ejecta. To determine 
the cause of the steepening of the BVR bands, further analysis of 
the photometry and spectroscopy was carried out and is detailed in 
Section|5] 



3.2 Colour Curves 



Figure 2. UBVRIJHK light curves of SN 2004et, from soon after core- 
collapse to +414 days post explosion. The light curves are constructed using 
the data in Table[T]and Tableland have not been corrected for reddening. 



The extinction 
to be 0.41 ± 



to SN 2004et, 
0.07 ma g us ing 



E(B-V) was estimated 
Na ID equivalent widths 



dZ witter. Munari & Morettil 20041) . Reddening corrected colour 
curves for SN 2004et using the optical and NIR data detailed above 
are shown in Figure [4] and Figure [5] The optical colour curves are 
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Table 3. Table of properties of our sample of type IIP SNe. 



Figure 3. BVR light curves of SN 2004et during the tail phase. The solid 
lines indicate the decline rates calculated from ~ 1 36-300 days post explo- 
sion, while the dashed lines indicate the decline rates calculated using the 
data from ~ 295^115 d. A steepening of the BVR light curves is seen while 
the slope of the / band light curve becomes marginally less steep. 



compa red to those o f other IIP SNe : SN 1999em (Leonard et al. 
l2002ah . SN 2003hn dOlivare j 120081; IKrisciunas etall |2009h. SN 
2005cs jPastorello et al.l2009h and SN 1990E dSchmidt et all 19931 ; 
iBenetti et alj 1994) . The explosion epoch, distance and extinction 
for each SN is detailed in Table [3] Few data are available in the 
NIR for type IIP SNe, particularly at late times so in this case 
for comparison we can only use SN 1999em, SN 1987 A, SN 
200 5cs and SN 2002h h. SN 2002hh was a highly reddened SN 
and IPozzq et"al] ( 2006) adopted a two component model for the 
extinction with total A v of 5 mag, which we use in Figure|3] 

As is expected, SN 2004et became monotonically redder af- 
ter explosion for ~ 100 days as it cooled. This can be seen in both 
the optical and NIR comparison colour curves. Later on when SN 
2004et enters the nebular phase, the colour curves become bluer. 
The colour evolution of SN 2004et is most similar to 'normal' type 
IIP SNe, SN 1999em and SN 1990E, while SN 2005cs showed 
a red peak at ~ 100 days, characteristic of low-luminosity SNe 
dPastorello etai]|2004 [2009^ . The V-H, V-K and J-K colours for 
SN 2004et are more blue at ~ 300 days than for SN 1987 A and SN 
2002hh. The V-J colour evolution appears to be similar to that of 
SN 1987A and SN 2002hh, which suggests the difference lies in the 
H and K bands being more blue than for the other SNe. This could 
be due to weaker features of the Brackett series of H in the spec- 
tra of SN 2004et, but we do not possess spectra of a high enough 
signal-to-noise to quantitatively determine if this is the cause. 

3.3 Bolometric light curve and 56 Ni Mass 

The bolometric light curve of SN 2004et was computed by cor- 
recting the observed magnitudes for reddening and converting to 
flux values before integrating the combined flux from the UBVRI- 
JHK bands exclud ing the overlapping regions of the filters (see also 
IValenti et all2008h . This integrated flux could then be converted to 
a luminosity using a distance of 5.9 Mpc (see Section[T} and the ex- 
tinction detailed in the previous section. Our UBVR1JHK data are 
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1987A 


2446850.5 


49.9±1.2 x 10~ 3 


0.19 
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1990E 


2447932.0 


18.2±2.4 


0.48 


2,3 


1999br 


2451278 


14.1±2.6 


0.02 


4 


1999em 


2451475.6 


11.7±1.0 


0.06 


5 


1999gi 


2451518.3 


10.0±0.8 


0.21 


6 


2002hh 


2452577.5 


5.9±0.4 


see text 


7, 8 


2003hn 


2452857.0 


17.8±1.0 


0.19 


9, 10 


2004A 


2453011.5 


20.3±3.4 


0.06 


11, 12 


2004et 


2453270.5 


5.9±0.4 


0.41 


13 


2005cs 


2453549.0 


7.1±1.2 


0.05 


14 


2006my 


2453943 


22.3±2.6 


0.03 


15 
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Figure 4. Optical colour evolution of SN 2004et compared with those of 
other type IIP SNe; SN 1999em, SN 2003hn, SN 2005cs, SN 1990E and 
SN 1987A. 



supplemented by optical data from lMisra et alj J2007I) . The bolo- 
metric luminosity was only calculated for epochs where V band 
data were obtained. If data were not obtained in the other bands 
at the same epoch, the contribution was calculated by interpolating 
the data from adjacent nights. 

Figure [6] shows a comparison with other type IIP SNe, some 
of which also include contributions from NIR data. Two pseudo- 
bolometric light curves of SN 2004et are shown, one that includes 
the contribution from the NIR bands (filled squares) and one that 
does not (open circles). When the NIR is included the luminosities 
of the plateau and the tail phase are significantly increased. The 
bolometric light curve of SN 2004et is one of few type IIP SNe 
that includes the contribution of NIR light, which is seen to make 
a significant contribution (up to 50 per cent). The percentage flux 
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Figure 5. Near infrared colour evolution of SN 2004et compared with those 
of SN 1999em, SN 2005cs and the peculiar SN 1987A. 



Figure 7. Flux contribution from NIR bands as a percentage of the total flux 
from optical and NIR bands for a selection of type IIP SNe with NIR data. 
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Figure 6. Comparison of bolometric light curve of SN 2()04et with other 
type IIP SNe. 



contribution of the JHK bands compared to the total flux from the 
optical to the NIR bands for SN 2004et and three other type IIP 
SNe, SN 1987 A, SN 1999em and SN 2005cs is shown in Figure|7] 
The contribution of NIR flux to the total flux of SN 2004et is most 
similar to the 'normal' type IIP, SN 1999em. It can be seen that 
particularly at early times the NIR flux contribution of SN 1987 A is 
significantly greater than that of other type IIP SNe due to the faster 
cooling of the ejecta of this SN. When no NIR data is available for 
a type IIP SN, the NIR contribution of SN 1987A is sometimes 
used as a correction to t he bolometric ligh t curve of SNe due to its 
excellent coverage (e.g. iMisra et alj[2007l) . Figure [7] demonstrates 
however that using a correction based on SN 1987A could lead 



to an overestimate of the bolometric luminosity, especially at very 
earl y times. 

ISahuetaI1 ( l2006l) computed a UBVRI bolometric light curve 
of SN 2004et and found it to have one of the highest luminosi- 
ties among t ype IIP SNe . They did not apply any correction for 
the NIR flux. lMisraet al. (2007) constructed a UVOIR bolometric 
light curve using a correction for the lack of NIR data, which was 
obtained by a comparison with the bolometric light curve of SN 
1987A. They used a smaller distance estimate to SN 2004et ( 5.5 
±1.0 Mpc), which results in smaller luminosity values. Despite 
the smaller distance used, their bolometric light curve displayed 
higher luminosities at early times than the light curve of SN 2004et 
shown in Figure [6] The peak luminosity o f the UBVR IJHK bolo- 
metric light curve of SN 2004et given in IMisra et al.l J2007t) was 
log L ~ 42.55 erg s compared to 42.44 erg s -1 for our data in 
Figure [6] This difference at early times is most likely due to their 
overestimate of the NIR luminosity from the comparison with SN 
1987A. 

iBersten & Hamuvl |2009) calculated the bolometric correction 
(BC) to transform from V band magnitudes of SN 1987A, SN 
1999em and SN 2003hn to UBVRIJHK bolometric luminosities 
and provided parametrised corrections as a function of SN colour 
for the average of these three SNe. Here we present the BC for four 
type IIP SNe with extensive optical and NIR coverage as a func- 
tion of time. These SNe are SN 1987A (a peculiar type IIP), SN 
2005cs (a low luminosity SN), and SN 1999em and SN 2004et, 
which both show a very simi lar bolometric light curve shape. We 
have rearranged equation 1 of lHamuvl d2003l) to produce this equa- 
tion for the BC at any given time t, 

B(t) = -2.5 log w L(t) - Vit) + A tot jy) + 5 log l0 D - 8.14 (1) 

where L(t) is the bolometric luminosity in ergs s , V(t) is the 
V band magnitude, A, om / is the total V band extinction toward the 
SN, D is the distance to the SN in cm and the constant is to convert 
from Vega magnitudes to values of luminosity in cgs units. Using 
the bolometric light curve for each SN, which has been constructed 
using the method detailed above, the BC as a function of time can 
be parametrised by fitting a polynomial to the data. The polyno- 
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Figure 8. Bolometric correction from the V band (upper panel) and R band 
(lower panel) to the UBVRIJHK bolometric light curve as a function of 
time for SN 1987A, SN 1999em, SN 2004et and SN 2005cs during the 
photospheric phase. The bolometric correction of each SN is fitted with a 
lower order polynomial (solid lines). 
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Figure 9. Bolometric correction from V band (upper panel) and R band 
(lower panel) to the UBVRIJHK bolometric light curve as a function of 
time for SN 1987A and SN 1999em during the early nebular phase. The 
bolometric correction of each SNe is fitted with a linear fit (solid lines). 



mial coefficients for the V and R bands during the plateau phase are 
given in Table lAll and Table I A2l respectivelv, 

The BC relative to the V and R bands during the photospheric 
phase are shown in Figure[8] For the V band, the scatter of the BC 
for each SN is relatively large with SN 1987A showing a signifi- 
cantly different evolution during the plateau phase. For the R band 
however, the evolution of the BC is much more homogeneous for 
the SNe, particularly after ~ 50 days and we suggest that the R band 
BC may be more suitable for computing the luminosities of a SN 
when there is ambiguity about the subtype of the type IIP SN being 
studied. 

The coefficients a, of the polynomial for each SN and band 
can be used in Equation [2] for each epoch of data to give the BC as 
a function of time, 

n 

BC(t) = Y j a,f (2) 

where n is the number of coefficients for the particular band 
and SN being studied. Equation[3]can then be used to give the bolo- 
metric luminosity of a SN with only single band photometric data, 

logwUt) = -0.4 [1/(0 - A ma ,(V) + BC(t) - 5 log l0 D + 8.14] (3) 

where the parameters have the same meaning as detailed 
above. Equation [3] is independent of the choice of zero point for 
the Vega magnitude system and is dependent only on the values of 
the magnitude, distance, extinction and time since explosion. 

Figure [9] shows the BC for the V and R bands to the UB- 
VRIJHK bolometric light curves of SN 1987A and SN 1999em 
during the early nebular phase (~ 120-200 days). SN 2004et and 
SN 2005cs did not have sufficient data at NIR wavelengths at 
these epochs to be included in the plot. The BC is relatively 
similar for both SN e, which is consistent with the results of 
lHamuv et all ( 1200 ll) and iBersten & Hamuvl d2009h . By applying 
the parametrised BC given in Table [A3] to data of other SNe with 
sparser photometric coverage during the early nebular phase, an es- 



timate of the UBVRIJHK luminosity can be calculated and hence 
an estimate of the ejected 56 Ni mass of a less well studied type IIP 
SN can be determined. We find an average BC during early nebular 
phase of 0.33 + 0.06 mag using SN 1987A and SN 1999em, which 
is systematically larger than the value of 0.26 ± 0.06 mag obtained 
by lHamuv et al.l |2001). The probable cause of this discrepancy is 
the different methods used to compute the bolometric light curve; 
we integrated the c ombined flux over th e UBVRIJHK bands as de- 
tailed above, while lHamuv et al" liooil) fitted Planck functions at 
each epoch to the reddening-corrected BVIJHK magnitudes. 

Information on the bolometric light curve is useful for con- 
straining the amount of radioactive material ejected in the explo- 
sion. We estimated the 56 Ni mass of SN 2004et using three meth- 
ods; a comparison w ith the luminosity of SN 1987A, the steepness 
parameter method oflElmhamdi et al .1 (12003) and the tail luminosity 
method of lHamuvl d2003l) . 

The 56 Ni mass can be estimated by a comparison of the bolo- 
metric luminosities of SN 2004et and SN 1987 A during the early 
nebular phase, before the possible formation of dust (~ 120-250 
days post explosion). T he 5S Ni mass of SN 1987A was 0.075 ± 
0.005 M Q dArnettlll996l) and based on a comparison of their UB- 
VRIJHK luminosities, we find SN 200 4et to have a mas s of 56 Ni 
of 0.057 ± 0.03 M G . Using this method lSahu et all d2006l) found a 
56 Ni mass of 0.048 ± 0.01 M Q . Given that for SN 2004et they did 
not include the NIR contribution, their value is a lower limit for the 
56 Ni mass. 

Despite the lack of good data coverage at the transition 
from the plateau to the tail phase, we used the equations of 
lElmhamdi. Chu gai & Danziged 12003) that correlate the steepness 
of the rate of the decline from plateau to tail phase of the bolomet- 
ric light curve and the mass of 56 Ni. The smaller the 56 Ni mass pro- 
duced in the explosion, the steeper the decline rate is. Using this 
method, we estimated the steepness parameter to have a value of 
0.068, which gives a 56 Ni mass of 0.057 ± 0.02 ISahu et al. 
(2006) found a 56 Ni mass of 0.062 ± 0.02 M Q and iMisra et~ 
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(2007) a 56 Ni mass of 0.056 ± 0.016 M using this steepness 
method. 

The 56 Ni mass of SN 2004et can also be estimated using the 
assumption that the tail luminosity of SN 2004et is dominated by 
the radioactive decay of 56 Ni to 56 Co (kamuy 2003) an d is found to 
have a value of 0.05 ± 0.02 M Q using the BC of lHamuv et al.l J200lh 
of 0.26 ± 0.06 mag. When our BC, calculated during the early neb- 
ular phase of 0.33 ± 0.06 mag is used, we obtain a 56 Ni mass of 
0.06 ± 0.02 M . The final assumed 56 Ni mass for SN 2004et used 
in this paper is taken as the average of the three methods listed 
above, which gives a value of 0.056 ± 0.04 M Q . 



4 SPECTROSCOPY 
4.1 Optical spectroscopy 

Optical spectroscopic observations of SN 2004et were carried out 
at the Mt. Ekar 1.82-m Copernico Telescope (Asiago, Italy) with 
AFOSC and at the TNG using DOLORES, as detailed in Table 
[4] The data were trimmed, bias corrected, and a normalised flat- 
field was applied before the spectra were extracted using tasks in 
iraf's ctioslit package. Wavelength calibration was performed us- 
ing arc lamp spectra. The instrumental response curves as a func- 
tion of wavelength were determined using spectrophotometric stan- 
dard stars. The spectra were then flux calibrated to the photometry 
taken at the closest epoch. The resolution of each spectra was esti- 
mated by measuring the full width half maximum (FWHM) of the 
night sky lines and the values are given in Table [4] 

An optical spectrum of SN 2004et taken at 56 days post explo- 
sion is shown in Figure [10] The key features of the plateau phase 
of the SN are marked. The identification of the lines was made by 
comparing t he spectrum to the tem plate plateau ph ase spectrum of 
SN 19 99em dLeonard et alj2002al) and SN 1998 A jPastorello et~ai1 
l2005h . for which detailed line identifications were carried out. The 
P-Cygni profile that is characteristic of a fast expanding ejecta is 
most clearly visible in Ha, but can also be seen in other elements 
marked in the spectrum. 

Figure QTJ shows the optical spectroscopic evolution of SN 
2004et during the photospheric p hase, from +24 to + 1 15 days. The 
spectra at +24 d and + 31 d from lSahuet all (2006) were retrieved 
from the SUSPEC10 archive to show the early time evolution of 
SN 2004et. As the photospheric phase progresses, the blue region 
of the spectra becomes dominated by metal lines, which add to the 
increasing absorption troughs of the H B aimer series. The absorp- 
tion troughs of the H lines become deeper as the photospheric phase 
progresses before decreasing as the SN enters the nebular phase. 

An optical nebular phase spectrum, taken at +401 days post 
explosion is shown in Figure [12] The most prominent emis- 
sion lines have been marked following the line identifications of 
iLeonard et al. (2002a) for SN 1999em. The spectrum is seen to be 
dominated by forbidden emission lines such as [O i] 6300, 6364 A, 
[Fen] 7155 A, [Can] 7291, 7324 A and a strong Ho- emission line. 
The Ho- and [Can] 7291, 7324 A doublet have comparable lumi- 
nosities. The individual components of the [Can] doublet are not 
resolved. Figure [T3l shows the overall spectroscopic evolution of 
SN 2004et d uring the nebular phase. The spectrum of +464 days 
is taken from lSahu et al.l d2006l) . At early nebular phases, it is seen 
that the most prominent feature is the Ho- emission line. The blue 
part of the nebular spectra are dominated by Fe n and [Fe n] lines. 




Figure 10. Optical spectrum of SN 2004et at 56 days post explosion (JD 
53326.0). The spectrum has been corrected for the recessional velocity of 
the host galaxy of 48 km s~ 1 . The prominent features are labelled. 



http://bruford.nhn.ou.edu/~suspect/index 1 .html 



lUtrobin & Chugall d2009» investigated possible asymmetry in the 
structure of SN 2004et by fitting three Gaussian components to 
each of the lines of the [Oi] 6300, 6364 A doublet in a +301 
day spectrum. Their results showed a non-negligible contribution 
from asymmetric components and from this, they inferred a bipolar 
structure for SN 2004et, whi ch has also bee n shown for SN 1987 A, 
SN 1999em and SN 2004dj dUtrobinll2007l) . 

A comparison between an optical photospheric spectrum of 
SN 2004et and other IIP SNe, SN 1999gi, SN 1999em, SN 2004A 
and SN 2005cs, all taken at ~ 40 days post explosion i s shown in 
Figure [141 The spectrum of SN 2004et is taken from ISahu et al.l 
(2006), while for SN 2004A see Appendix B. The other spectra 
are from the references given in Tabled SN 2004et has the most 
blue-shifted absorption minima of all the SNe of this sample, which 
suggests a higher expansion velocity for the outer ejecta. However 
it will be seen in Section [43] that the photospheric expansion ve- 
locity of SN 2004et calculated from weak metal lines is similar to 
that of SN 1999em. 



4.2 Near infrared spectroscopy 

NIR spectroscopic observations were obtained at 22 epochs from 
+7-65 days post explosion at the 1.08-m AZT-24 telescope Campo 
Imperatore Observatory (Italy) with SWIRCAM, as detailed in Ta- 
ble[5] A late time spectrum was obtained with the TNG+NICS on 
25 July 2005 (+306 days post explosion). 

The reduction of the NIR spectra was carried out using iraf's 
onedspec package. Pairs of spectra (labelled spectrum A, spectrum 
B) were obtained immediately after each other with the SN at two 
different positions along the slit. Each pair of images were sub- 
tracted from each other (A-B) so that the sky background was re- 
moved from the 2-D spectrum of the SN. Three pairs of images 
were obtained for each epoch and for each grism. These subtracted 
pairs were then added together to produce a 2-D spectrum of higher 
signal-to-noise. From these combined images, the SN spectrum 
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Figure 11. Optical spectroscopic evolution of SN 2004et during the photo- 
spheric phase with days marked since the explosion epoch (JD 2453270.5 
). The spectra have been corrected for the recessional velocity of the host 
galaxy of 48 km s . Some key emission and absorp tion features are la- 
belled. The spectra of +24 d and +31 d are taken from lSahu et alj 12006). 

was then extracted and any residual background that remained in 
the wings of the aperture was also subtracted off. 

No arc lamp spectra were available so the sky lines in the orig- 
inal non-subtracted images were used to perform the wavelength 
calibration for the early time spectra. This was carried out by ex- 
tracting along the same line as the SN but in the other image of the 
pair so that only the background along that line remained. These 
sky spectra were then used as 'arcs' due to the well known features 
of atmospheric OUT ions. A further check of the wavelength cali- 
bration was carried out using the telluric features present in the SN 
spectra. 

The late time spectrum from the TNG was obtained using the 
low resolution Amici prism with a wavelength range of 0.8-2.5 fim. 
The resolution of this spectrum is not given in Table[5]because the 
Amici prism has a nearly constant resolving power and so the dis- 
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Figure 12. Optical spectrum of SN 2004et at 401 days post explosion. The 
spectrum has been corrected for the recessional velocity of the host galaxy 
of 48 km s . The prominent features are labelled. 

persion varies by more than a factor of 3 over the spectral range. 
Due to this very low resolution the arc lamp lines are blended so 
a look-up tablqj detailing pixel coordinate against wavelength is 
used to perform the wavelength calibration. The pixel coordinates 
are converted to wavelengths using this table and then a further cor- 
rection is performed to the fit using prominent emission lines of the 
SN. The SN spectrum was divided by a G-type telluric standard 
star to remove the strong NIR telluric features. The instrumental 
response of the system was found by comparing the standard star 
observation to that of a solar type spectrum. Then the SN spec- 
trum was multiplied by this response function to obtain the flux 
calibrated spectrum. Finally the flux calibration of the SN spectra 
was checked against the NIR photometry from the same epoch and 
adjusted if necessary. 

3 http://www.tng.iac.es/instruments/nics/spectroscopy.html 
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Figure 13. Optical spectroscopic evolution of SN 2004et during the nebular 
phase with days since the explosion epoch marked next to spectra. Emission 
lines are seen to be the dominant features during this phase. 
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The spectral evolution of the NIR HK band (-1.45-2.2 yum) 
is shown in Figure [16] Five members of the Brackett series of H 
have features that are visible in the HK band spectra. The HK band 
does not show any major evolution between +7 days and +63 days. 
Like the IJ band spectra, the resolution of the spectra is low so not 
many spectral features other than the H Brackett series, Hei and 
possibly Mg i can be identified. The feature ide ntification of both 
the IJ and HK bands was carried out following lFassia e t al. (2001) 
and lGerardv et al.l J2000h . 

Figure [T7] shows a comparison of two photospheric phase I J 
band spectra of SN 2004et with other type IIP SNe for which 
NIR spectra were a vailable at roug hly coeval epochs, S N 1997D 
( iBenetti et all 1200 ll). SN 199 9em lHamuv et all 1200 ll) and SN 
2005cs jPastorelloetalJl2009h . To ensure a realistic comparison 
between spectra of different resolutions, the spectra have been 
smoothed to the resolution of the worst spectrum. The lowest res- 
olution spectra are those of SN 2004et with a resolution of 38 A. 
The comparison spectra have been convolved with a Gaussian with 
a FWHM equal to this resolution and then the smoothed spectra 
have rebinned to a common dispersion relation. Note that the tel- 
luric features in the spectrum of SN 2005cs have not been well 
removed because the telluric standard spectrum w as taken a few 
days after that of SN 2005cs dPastorello et al]2009h . 

At ~ 20 days, the spectra of SN 2004et and SN 1999em show 
mainly the features of the Paschen series of H. For the later epoch, 
we note that SN 1997D and SN 2005cs had lower expansion veloc- 
ities and so the separate components making up the feature at ~ 1.1 
/im are identifiable, unlike in the spectra of SN 2004et. The NIR 
spectra of this selection of type IIP SNe do show many similarities, 
with the H Paschen series being particularly prominent in all the 
spectra. Four combined optical and NIR spectra of SN 2004et are 
shown in Figure [T8l Some of the o ptical spectra used in this figure 
were taken from Sahu et al.l d2006h due to their larger wavelength 
coverage and matching epochs with our NIR data. 

The late t ime NIR spectrum of S N 2004et is com pared to that 
of SN 2005cs dPastorello et~ai1l2009l) and SN 2002hh dPozzo et al.l 
120061) in Figure[l9l Th e line identifications of lHamuv et alj feOOlh 
and lpozzo et al. (2006) were adopted for the nebular spectrum. As 
in the photospheric phase spectra of SN 2004et the prominent fea- 
tures are those of the Paschen series of H, along with the other 
heavier elements marked in Figure[T9] Like in the spectrum of SN 
2002hh, the first overtone band of CO is also seen in SN 2004et 
at wavelengths longer than 2.3 fjm. This is known as a indicator of 
dust formation and is discussed in more detail in Section[5] 



Figure 14. Comparison of coeval optical photospheric spectra of type IIP 
SNe at ~ 39 days post explosion. The spectrum of SN 2004et is taken from 
ISahu et alj feOOtj) and the spectrum of SN 2004A is from Appendix B. The 
other spectra are taken from the references given in Tabled 



Figure[T5]shows the spectral evolution of the NIR IJ (~ 0.85- 
1.35 (im) band during the photospheric phase. The spectral features 
show the characteristic P-Cygni profiles of a fast moving ejecta, 
which are also visible in the optical spectra during the photospheric 
phase. Some lines of the Paschen series of H are seen, along with 
Ci, He 1 and Oi lines. The wavelength region between ~ 1.8 yum 
to just before 2 fim was lost due to strong atmospheric absorption 
in this region. The spectra are low resolution so not many spectral 
features can be determined unambiguously. The resolution of each 
spectra was calculated from the slit width and the dispersion and 
these values are quoted in Table[5] 



4.3 Expansion velocity 

The expansion velocities of Hor (6563 A), H/? (4861 A), Sc 11 (6246 
A) and the average of the Fen triplet at 4924, 5018, 5169 A are 
shown in Figure l20l as determined by fitting a Gaussian to the ab- 
sorption profile of each line and measuring the position of the min- 
imum. As is typical for type IIP SNe, Ho- gives higher velocities, 
indicating line formation at larger radii. Lines with lower optical 
depths such as Fen 4924, 5018, 5169 A are better indicators of the 
photospheric velocity. Even better in this respect is the Sc n 6246 A 
line, which has an even lower optical depth and hence is the most 
suitable line for determining the expansion velocity of the photo- 
sphere. 

Figure|2T]shows the photospheric velocity estimated from the 
Sc 11 (6246 A) line of SN 2004et comp ared with the velocit ies for 
three oth er type IIP SNe, SN 1999br dPastorello et~aT1l2004l). SN 
1999em ^Leonard et all 12002 j) and SN 2005cs dPastorello et"ai1 
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Table 5. Log of near infrared spectroscopic observations. 



Date 


JD (2450000+) 


Phase* (days) 


Telescope +instrument 


Grism 


Range (pm) 


Resolution (A) 


29/09/2004 


53277.9 


7.4 


A7T-?4+SWlRrAM 


/ + / H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38 72 


10/09/2004 


53278 8 


8.3 


A7T-24+SW1RCAM 


/ + / H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38 72 


01 /I 0/7004 


53279.9 


9.4 


A7T-24+SWlRrAM 


l + j H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38, 72 


09/1 0/9004. 


JJiOv/.O 


10 3 


A 7 T- 9 4. + S W TR C A M 


1+ J H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38 72 


0V1 0/9004 


53281 9 


114 


A 7 T- 9 4. -4- S W TR r A M 


1+ J H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38 72 


04./ 1 0/9004. 




12 3 


A 7 T- 9 4. -4- *n W TR C A M 


1+ J H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38 72 


0S/1 0/9004. 


JjZOJ.7 


13 4 


A7T-94.+*NWTRr AM 


1+ J H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38 72 


07/1 0/9004. 


JJiOJ .o 


15 3 


A7T-94.+SWTRT AM 


1+ J H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38 72 


0S/1 0/9004. 


JJZ.OU.O 


16.3 


A 7 T- 9 4. 4- S W TR C A M 


1+ J H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38 72 


09/10/2004 


53287.8 


17.3 


AZT-24+SWIRCAM 


I+J, H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38,72 


18/10/2004 


53296.8 


26.3 


AZT-24+SWIRCAM 


I+J 




0.84-1.32 




38 


20/10/2004 


53298.9 


28.4 


AZT-24+SWIRCAM 


I+J, H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38,72 


21/10/2004 


53299.8 


29.3 


AZT-24+SWIRCAM 


I+J 




0.84-1.32 




38 


22/10/2004 


53300.8 


30.3 


AZT-24+SWIRCAM 


I+J 




0.84-1.32 




38 


23/10/2004 


53301.7 


31.2 


AZT-24+SWIRCAM 


I+J, H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38,72 


25/10/2004 


53303.8 


33.3 


AZT-24+SWIRCAM 


I+J, H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38,72 


27/10/2004 


53305.8 


35.3 


AZT-24+SWIRCAM 


I+J, H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38, 72 


17/1 1/2004 


53326.8 


56.3 


AZT-24+SWIRCAM 


I+J, H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38,72 


21/11/2004 


53330.8 


60.3 


AZT-24+SWIRCAM 


I+J, H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38,72 


24/11/2004 


53333.7 


63.2 


AZT-24+SWIRCAM 


I+J, H+K 


0.84- 


-1.32, 1.45- 


-2.38 


38,72 


06/12/2004 


53345.7 


75.2 


AZT-24+SWIRCAM 


I+J 




0.84-1.32 




19 


25/07/2005 


53576.6 


306.1 


TNG+NICS 


AMICI prism 




0.8-2.5 




(see text) 



* relative to the epoch of date of explosion (JD = 2453270.5) 




Wavelength (A) Wavelength (A) 



Figure 15. IJ band spectral evolution of SN 2004et during the photospheric 
phase from +8 days to +75 days post explosion. The positions of the most 
prominent spectral features are marked. 



Figure 16. HK band spectroscopy during the photospheric phase from +7 
days to +63 days post explosion. The positions of the most prominent spec- 
tral features are marked. 



2009). The Sen velocities of the comparison objects were taken 
from the literature. The photospheric velocity of SN 2004et is seen 
to have a similar evolution to that of SN 1999em. The velocities of 
SN 2004et and SN 1999em are, as expected, higher than the veloc- 
ities of the low luminosity type IIP SNe 1999br and 2005cs. 



5 DUST FORMATION 

Dust can form in SN ejecta as they cool during the neb- 
ular phase. One of the first signs of dust formation is the 
formation of molecules containing carbo n and silicon, which 
are the precursors to dust formation fcjerardv et al.l l200Ct 
Spyromilio, Leibundaut & Gilmozzi 2001). The first overtone of 
CO emission is expected at NIR wavelengths longer than ~ 2.3 /jm 
before the onset of dust. This signature of dust formation is seen 
in the spectrum of SN 2004et taken at ~ 300 days post explosion, 
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Figure 17. Comparison of the IJ band spectra of SN 1999em, SN 1997D 
and SN 2005cs with two spectra of SN 2004et at epochs roughly matching 
the comparison SNe. The positions of the most prominent spectral features 
are marked. 



Figure 19. The late time nebular (~ 300 days) N IR spectrum of SN 20 04et 
is compared to that of SN 2005cs at ~ 281 da ys IPastorello et alj2009l) and 
SN 2()02hh at + 266 days jPozzo et alj200rj) . 
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Figure 18. Combined optical and NIR spectra of SN 2004et. For the epochs 
at +55 days and +30 days, the NIR spectra were obtained one day later 
than the optical observation. The NIR spectra were combined with optical 
spectra taken from Sahu et al. (2006) and this paper. 



which is shown in Figure [19] The SN 2002hh also shows a simi- 
lar CO emission at these wavelengths, unlike SN 2005cs where no 
significant CO emission band is seen. No earlier NIR nebular spec- 
tra of SN 2004et are available to constrain the epoch of molecular 
formation and no attempt is made to estimate the quantity of CO 
formed due to the low resolution of the spectrum. The presence of 
the CO molecular band has been seen in a good fraction of type 
IIP SNe with observati ons at NIR wave lengths such as SN 1987 A 
jSDvromilioetal.lll988h, SN 1995ad JSpyromilio & Leibundgul 
Il99fj) . SN 1998dl dSpyromilio. Leibundgut & Girmozzill200ll) . SN 
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Figure 20. Velocity evolution of some prominent lines of SN 2004et, mea- 
sured from the minima of their absorption profiles. 



1999em jSpyromilio, Leibundgut & Gilmozzill20oih , SN 2002hh 
dPozzo et all2006h and now SN 2004et. This suggests that type IIP 
SNe must cool significantly within a few hundred days of explosion 
to allow first molecular formation and then dust formation. 

Other indicators of dust formation include a blueshift of the 
spectral emission lines of elements such as Hcf and [On] 6300, 
6364 A, an optical luminosity decrease, a NIR excess at late times 
or thermal emission from newly formed dust grains at MIR wave- 
lengths. 

A blueshift of the peaks of emission lines is caused by resid- 
ual opacity in the central part of the ejecta due to dust particles 
and results in an attenuation of the red wing of the line profile. 
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Figure 22. Evolution of the line profiles of Her for SN 2004et (left panel) 
and for SN 1999em (right panel). The vertical line in the panels corresponds 
to the zero velocity of Ho- at 6563 A. 
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Figure 21. Comparison of the expansion velocities of SN 2004et and three 
other type IIP SNe: SN 1999br, SN 1999em and SN 2005cs, using the Sen 
(6246 A) line. The Sc n line was not visible in the three earliest spectra of 
SN 2004et and so the expansion velocity for these epochs has been esti- 
mated from 0.95 times the velocity of the average of the Fe n 4924, 5018, 
5169 A lines. 

ISahu et alj fim& noted that this blueshift in the peaks of the Ha 
and [ Oi] 6300, 6363 A occurred after ~ 300 days. Using the spec- 
tra of lSahu et al.1 12006) and the additional spectra published in this 
work, the blueshift of the emission lines was quantified. The size of 
the shift can be determined by measuring the wavelength position 
of the centroid of the Ho- emission line. For the spectra between 
163 days and 300 days post explosion, the peak of the Ha emission 
line was found to be at +280 ± 50 km s . The spectrum at 313 
days and all subsequent spectra out to the last at 464 days show a 
constant blueshift to -137 km s -1 . 

Figure [22] and Figure [23] show the temporal evolution of the 
Ha and [Oi] 6300, 6364 A line profiles in SN 2 004et compared 
with those of SN 1999em. lElrnhamdi et al.1 fc003h showed for SN 
1999em that there was an observable blueshift of the emission lines 
at ~ 500 days, which they suggest is caused by dust formation. We 
measured the rest frame position of the Ha peak in SN 1999em to 
be +182 km s -1 at 466 days shifting t o -91 km s -1 at 642 days using 
spectra taken from lElmhamdi et alj |2003). Dan ziger et alj J 1 99 lh 
showed that for SN 1987 A, an observable blueshift was seen in 
the emission peaks of the [0 1] doublet, beyond approximately 400 
days post explosion. This can again be measured from the Ha peak 
as +164 km s" 1 at +440 days, -156 km s -1 at +501 days and then 
increasing to -320 km s -1 at +628 days. The larger blueshift in 
the emission lines of SN 1987A compared to SN 1999em can be 
attributed to either a larger mass of dust forming in the ejecta or to 
a geometrical effect. The emission lines of SN 2004et experienced a 
blueshift at an earlier epoch than either SN 1987A and SN 1999em, 
and had intermediate values to those SN 1987A and SN 1999em. 

A decrease in the luminosity of the optical bands at late time is 
also thought to be indicative of scattering of optical photons off dust 
particles formed in the inner or outer envelope. Figure [3] showed a 
steepening in the slope of the optical BVR bands, as detailed in 
Section [3~T| and suggests that the steepening of the slope occurs at 
~ 300 days. This epoch of 300 days for SN 2004et is consistent 




Figure 23. Evolution of the line profiles of [0 1] 6300, 6364 A for SN 2004et 
(left panel) and for SN 1999em (right panel). The vertical lines in the pan- 
els correspond to the zero velocities of the [0 1] doublet at 6300, 6364 A 
respectively. 

with the dust formation rang e estimated from th e blueshift of neb- 
ular emission lines. Finally, iKotak et al.1 J2009h observed thermal 
emission from newly formed dust grains at MIR wavelengths and 
presented evidence for an IR echo from the interstellar dust of the 
host galaxy that manifests itself as a cold component of the spectral 
energy distribution of SN 2004et. 



6 ESTIMATE OF O MASS 

The mass of O in type IIP SNe can be estimated f rom an analy- 
sis of the [O i] 6300 6364 A lines dUomotol 1 19861 : Ichugail 1 1994 
iLi & McCravlll992l ; lElmhamdi et al.ll2003l) . The luminosity of the 
[Oi] doublet in SN 2004et can be compared to that of SN 1987A 
during the late nebular phase but before dust formation. Dust for- 
mation did not occur until at least 400 days post explosion for 
1987 A jDanzigeretalJl99l[) but for SN 2004et dust formation ap- 
peared to occur at ~ 300 days post explosion. To make a consistent 
comparison, spectra from ~ 285 days post explosion were used. 
The luminosities of the [0 1] 6300, 6364 A doublet along with the 
luminosities of what is thought to be the [Can] 7291, 7324 A dou- 
blet are given in Table[6]for a selection of type IIP SNe. The feature 
at 7300 A could also have a contribution from the [O n] 73 19, 7330 
A doublet, particularly as SN 2004et is likely to be an O-rich SN 
(see below). 

The s pectrum of SN 1987 A from 1987 December 09 was 
taken from lTerndrup et al.1 i 19881) . Spectra of SN 1999em and SN 
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2005cs were not available at the epoch of ~ 285 days post explosion 
so spectra before and after this epoch were used and an interpolated 
value was der i ved. Th e spectra of SN 1999em were taken from 
iLeonard et al.l |2002a) at +168 days and +313 d ays post explo- 
sion. T he spectra of SN 2005cs were obtained from Pastorel To et al.l 
d2009h at +281 days and +333 days post explosion. Flux calibration 
of the spectra used has been carried out by comparing to photome- 
try at the same epoch. 

The luminosity of the [Oi] doublet is lower for SN 2004et 
than for SN 1987A. The [Oi] doublet luminosity at late epochs 
is powered by y-rays being deposited in O-rich material and a 
relation between O m ass and line luminosity can be written as 
dElmhamdi et alj|2003h : 

1011 = 77 M ( } 

where M„ r is the 'excited' mass in which the y-ray deposition 
takes place, rj is the efficiency of the transformation of the energy 
deposited in the O mass into the luminosity of the [Oi] doublet, 
L[ i] is the O luminosity and L Co is the luminosity of 56 Co. The 
ratio of the mass of O in SN 2004et to SN 1987 A can be estimated, 
making the assumptions that 77 and the 'excited' mass are similar 
for both SNe. 

The Leo is directly related to the mass of 56 Ni, which was es- 
timated from the nebular phase bolometric light curve in Section 
[331 to be 0.056 ± 0.04 M . Th e O mass for SN 1987 A was found 
to be in the range, 1.2-1 5 M s dLi & McCravIll 9921 ; IChugaill 1994 
iKozma & FranssorJ 19981) . In SN 2004et, the [0 1] doublet luminos- 
ity is ~ 35 per cent lower than in SN 1987A and at the same time the 
56 Ni mass is ~ 25 per cent lower. The two differences almost match 
and following Equation [4] a similar O mass is derived. The values 
obtained for the O mass of SN 1999em and SN 2005cs are also 
given in Table|6] where the 56 Ni mass of SN 1999em was estimated 
to be 0.05 M from a comparison of its bolometric luminosity with 
that of SN 1987 A using t he Cepheid distance of 11.7 Mpc taken 
from lLeonard et alj d2003l). The 56 Ni mas s of SN 2005cs was taken 
as 0. 003 M fromlPastorello et al.l d2009l) . 

lUomotold 19861) estimated the minimum mass of O needed to 
produce the [0 1] emission lines using the equation: 

Af [0 i] = 10 8 F [Ol] »V- 28/r4 (5) 

where F[ ij is the [Oi] doublet flux in units of ergs s" 1 , D is 
the distance to the SN in units of Mpc and T4 is the temperature 
in units of 10 4 K. From Liu & DalgarnrJ d 19951) . the O temperature 
of SN 1987A at 300 days was ~ 4200 K. Assuming a similar O 
temperature for SN 2004et at a comparable epoch, the M[oi] of SN 
2004et was calculated for temperatures in the range 3500-4500 K. 
The minimum O masses determined using Equation [5] are given in 
Table |U For SN 1999em and SN 2005cs the O temperature was 
also assumed to be similar to that of SN 1987 A and the lower limit 
of the mass of O needed to produce the [0 1] emission lines are also 
estimated from Equation|5] 

The lower limit of t he O mass for SN 198 7 A i s consistent wit h 
the O mass range set bv lli & McCravl dl992l) and lChugail dl994l) . 
Equation [5] is most sensitive to the O temperature and this is the 
main source of uncertainty in the calculation. The mass of O in SN 
2004et is very similar to that of SN 1987A, while the O masses 
of the 'normal' type IIP, SN 1999em and the low luminosity SN 
2005cs are found using both methods to be lower. Despite the low 
luminosity of the [Oi] doublet seen for SN 2005cs, Equation [4] is 
inversely proportional to the L Co , which is also low for SN 2005cs 
and this results in similar O masses for SN 1999em and SN 2005cs. 
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Figure 24. Late time multi-band photometry of SN 2004et from +300 days 
post explosion to the latest Gemini North image taken on 2008 July 10 
(+1387 days). 

The similar O masses of SN 2004et and SN 1987A could suggest 
that they have similar main sequence progenitor masses, while SN 
2005cs and SN 1999em have lower O masses and could be asso- 
ciated with lower mass progenitor stars. IWooslev & Weaverl fl995) 
suggested that a main sequence star of mass of 20 M G would pro- 
duce an O mass of 1.5 M Q and so points toward relatively high mass 
stars as the progenitors of type IIP SNe. The discrepancies between 
direct imaging of SN progenitors and other estimates of the mass 
will be discussed in Section [8] for SN 2004et and more generally 
for type IIP SNe in Section|9] 



7 VERY LATE TIME PHOTOMETRY OF SN 2004et 

The late time optical and NIR HST data described in Section|2]and 
detailed in Table Q] and Table [2] are discussed here in more detail. 
The two epochs of HST data were obtained at +1019 and +1214 
days post explosion in five filters (two optical and three NIR). 
ICrockett et~al] d2009ri also obtained observations of SN 2004et with 
the WHT in four filters, BVRI on 2007 August 12 (+1054 days), 
along with a ground-based adaptive optics image with the Gemini 
North Telescope using Altair/NIRI on 2008 July 10 (+1387 days), 
which gave K * 20.0 mag. 

The late time evolution of SN 2004et at optical and NIR wave- 
lengths is shown in Figure|24] The HST magnitudes have not been 
converted to standard V1JHK filters because of the uncertainties in 
the transformation equations for SN spectra at such late phases. In 
this respect, the increase in the magnitude between the +1019 and 
+ 1054 days observations in the / band should not be regarded as 
significant. 

Ignoring this small increase, the trend between +1019 and 
+ 1214 days in all the bands is a levelling off of the magnitude de- 
crease. The light curves at early times are powered by the radioac- 
tive decay of 56 Co to 56 Fe but at around +1000 d ays other radioac- 
tive nu clides become important. In particular iFransson & Kozmal 
d2002h showed for SN 1987A that the decay of S7 Co became dom- 
inant after ~ 1100 days, while at epochs greater than 2000 days, 
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Table 6. Spectral information used in Equation U] and Equation[5] The spectra of SN 1987A and SN 2004et are both taken at ~ 285 days post explosion. The 
luminosities for SN 1999em and SN 2005cs were interpolated as detailed in Section[6] 



SN 


Luminosity [0 1] 6300, 6364 A 
( X 10 39 ergs-') 


Luminosity [Can] 7291, 7324 A 
(X 10 39 ergs-') 


Min. O mass from 
Equation|5](Mo) 


O mass from 
Equation|4](Mo) 


1987A 


2.0 ±0.1 


5.6 ±0.2 


0.3-1.1 


1.2-1.5^ 


1999em 


0.7 ± 0.4 


1.4 ±0.5 


0.1-0.4 


0.6-0." 


2004et 


1.6 ±0.4 


3.6 ±0.5 


0.2-0.9 


1.2-1.5 


2005cs 


0.05 ± 0.02 


0.15 ±0.02 


0.003-0.01 


0.8-0.9 



a Mass range estimate taken from lLi & McCravl jl992h and lChugail Jl994l) . 



44 Ti is the dominant source of energy. This could account for some 
of the flattening in the light curves between +1019 and +1214 days 
but not for the clear increase the magnitude that is seen in the Gem- 
ini A" band magnitude at +1387 days. 

The late time optical spectra (after ~ 800 days) of iKotak et al.l 

(2009) showed signs of an interaction between the ejecta and CSM 
in the wide boxy profiles of the Ha, [Can] and [Oi] emission lines. 
This may be the so urce for the additional luminosity at late times. 
IKotak et alj d2009h also suggested that there could be late time 
contributions to the luminosity from dust condensation in the cool 
dense shell left behind the reverse shock wave. The presence of in- 
teraction bet ween the SN ejecta and the progenitor CSM is further 
backed up by IKotak etHll ( 120091) by the discovery of a significant 
re-brightening in mid infrared Spitzer data. They also suggested a 
flux increase at optical wavelengths, which is however less con- 
vincing. They estimated the broad-band magnitudes of the BVRI 
bands by fitting a black body template to late time Keck spectra. 
Their magnitudes estimated from a spectrum at +1146 days are ~ 
0.5 mag brighter in the V band and ~ 0.3 mag brighter in the / band 
than the WHT images at + 1054 days. In addition, their data set does 
not include the HST photometry from +1214 days, which does not 
show any increase and are actually slightly fainter than those taken 
with HST at +1019 days. Therefore Figure l24l does not include 
the ma gnitudes estimated from the spectra presented in Kot ak et all 
d2009l) . 



8 PROGENITOR MASS: HIGH MASS OR LOW MASS? 

The focus of this work is to present and analyse the extensive data 
set of SN 2004et and therefore detailed modelling of this SN is out- 
side the purpose of this investigation. For most of the results on hy- 
drodynamic modelling reported in the following sections, we will 
refer t o work already published in the literature. lUtrobin & Chugail 
(2009) used a detailed 1-D hydrodynamical code to estimate the 
progenitor mass of SN 2004et. They modelled the bolometric light 
curve and the photospheric velocities of SN 2004et using hydrody- 
namical simulations to obtain an ejected mass of 24.5 ± 1 M Q and 
an estimate d progenitor mass on the main sequence in the range 
25-29 M^. IUtrobin&Chu gai (2009) suggested that the mass for 
SN 2004et obtained from the hydrodynamical modelling could be 
overestimated due to their use of a one-dimensional approximation. 
This approximation only artificially includes the effects of mixing 
between the He core and the H envelope and clearly cannot take 
into account any asphericity in the explosion. We stress that ev- 
idence for asymmetries have been fo und from polarimet ric mea- 
surements for five other IIP SNe (e.g. [Leonard et all 2006). 



lYound d2004l) performed a hydrodynamical parameter study 
of type II SNe and SN 2004et can be compared to the models in 
this study. SN 2004et is found to be most similar to model B, which 
is a 20 M main sequence mass star with a pre-explosion radius of 
431 R , explosion energy of 10 51 ergs and an ejected nickel mass 
of 0.07 M Q . The ejected mass taking into account mass loss and 
remnant mass is 16 M Q . The resulting observational parameters of 
this model are a plateau duration of ~ 100 days and an absolute 
V band magnitude of — 16.7 mag, which are both similar to the 
values measured observationally for SN 2004et. 

In Section [9] the ejected mass of SN 2004et is estima t ed us- 
ing the hydrodynamical equations of lLitvinova & Nadezhinl dl985l) 
and observational properties of SN 2004et during the plateau phase, 
to be 14 ± 6 M e . The mass of the neutro n star remnant of ~ 1.4 
M s dBogdanov. Rvbicki & Grindlay||2007l) and the mass loss due 
to stellar winds of ~ 1 M n , whi ch gives a main sequence mass of 



~ 16.4 ± 6 M e . Chugai 1 19941) used the nucleosynthesis models 
of IWooslev & Weaver! d 19951) to estimate the main sequence mass 
of SN 1987A to be ~ 20 M from their calculated O mass. The O 
mass range obtained for SN 2004et from the analysis of the nebu- 
lar spectra is seen to be of a comparable size to the O mass of SN 
1987A and hence we may guess a similar main sequence mass for 
SN 2004et of ~ 20 M G . 

IChevalier. Fransson & Nvmarkl J2006h showed that the radio 
and X-ray properties of SN 2004et were similar to those of other 
'normal' Type IIP SNe. SN 2004et had the highest radio luminosity 
of the sample studied. It was also shown using radio properties that 
the mass loss rate of ~ 2 x 10~ 6 M Q yr~' before explosion for SN 
2004et is consistent with the mass-loss rate expected for a progen- 
itor star of ~ 20 M Q . 

These estimates of the progenitor mass of SN 2004et suggest 
that the main sequence mass is toward the higher end of masses 
for Type IIP SNe. This seems consistent with the bolometric light 
curve, which has one of the highest luminosities of the sample of 
IIP SNe shown in Figure |6] The 56 Ni mass ejected by SN 2004et 
was found to have a value of 0.056 ± 0.04 M Q , which is relatively 
high for a type IIP SN. 

However dCrockett et aIT2 009) have studied pre-explosion im- 
ages of the site of SN 2004et and find a progenitor star with a 
main sequence mass of 8^j M Q . This value is lower than the es- 
timates determined from studying the explosion parameters. The 
discrepancies between the progenitor mass estimates for type IIP 
S Ne obtained u s ing th ese t wo approaches has already been noted 
by ISmartt et alj {2009) and lUtrobin & Chugail d2008t |2009|) . Sec- 
tion [9] discusses in more detail these discrepancies for a sample of 
ten well-studied type IIP SNe. 
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9 EXPLOSION PARAMETERS OF TYPE IIP SNE 

The progenitor and explosion properties of type IIP SNe can be 
studied in a number of ways. The observational properties of a 
SN such as the magnitude, expansion velocity and plateau length 
can be measured and models can then be used to determine the 
explosion parameters such as the ejected mass, explosion energ y 
and pre-explosion rad i us of the s tar (e.g. IZampieri et alj 120030 . 
iLitvinova & Nadezhinl ( 1983, 1985) used these observational prop- 
erties of SNe at mid-plateau as inputs t o hydrodynamic al equations 
to estimate their explosion parameters. lHamuvl Jjjolh determined 
th e physical properties o f 24 ty pe IIP SNe, using the equations 
of ILitvinova & Nadezhinl d 19851) but obtained much higher val- 
ues than those obtained by pre-explosion imaging, parti cularly for 
the eje cted mass with values in the range of 14-56 Mo-Eadezhin 
d2003D also derived the properties of a sample of ty pe IIP SNe, us- 
ing a subset of the observational data presented bv lHamuvl d2003l) 
and found values for the ejected mass of 10-30 M G , which are 
more consistent with those obtained using direct imaging of the 
pre-explosion stars. 

More complex hydrodynamical models have also been devel- 
oped for the analysis of the explosion parameters of type IIP SNe. 
SNe with good photometric and spectroscopic coverage have been 
studied such as SN 1987 A dUtrobmll 19931; iBlinnikov et aI.ll200Cl: 
Utrobml2004h , SN 1997D dChugai & Utrobinl200d:IZampieri et al] 
20031) , S N 1999br dZampieri et alj2003l;|Pastorello et al.l2004l). SN 



1999e m dUtrobinll2007l). SN 2003 Z ( Utrobi njChugai & Pastorellol 
|2007|). SN 2004et dUtrobin & Chugail 120091) and SN 2005cs 
dUtrobin & Chugaill2008l ; iPastorello et alj |2009). They have used 
models that include not only the explosion energy, pre-explosion 
radius and ejected mass as the input parameters but also the 56 Ni 
mass. This leads to a different set of relations for determining 
the physical parameter s of th e explosion than those de termined by 
ILitvinova & Nadezhml |l985). Zampieri 12005 1 |2007l) presented a 
systematic analysis of type IIP SNe using a 1-D, Lagrangian radia- 
tion hydrodynamics code and a semi-analytic code. The SN param- 
eters were estimated by performing a simultaneous comparison of 
the observed and simulated light curves, the evolution of line ve- 
locities and the continuum temperature. In some cases, these more 
complex models also find significantly higher masses than those 
obtained from direct imaging but for other SNe there is reasonably 
good agreement (cf Table[7jl. SN 1987 A shows excellent agreement 
bet ween the mass of its progenitor and hydrodynamical modell ing 
(e.g. lUtrobinlll993l ; lBlinnikov et al.ll200(]| - |Pastorello et alj2005h . 

ISmartt et al. ( 2009) presented estimates of the masses (and up- 
per mass limits) of the progenitors of 20 of the nearest IIP SNe in 
the last 10 years. These mass estimates came from deep high res- 
olution images of the SN field before explosion and gave an ini- 
tial mass from the progenitors in the range of 8-17 M Q . If it is 
assumed that a typical red supergiant loses around 1 M Q due to 
stellar winds and a ne utron star remnant is left behind with an av- 
erage mass of 1.4 M Q (Bogdanov, Rvbicki & Grindlav 2007), then 
the ejected masses would be in the range of 6-15 M Q , which are in 
general lower than those obtained from hydrodynamical modelling. 

Although the models of ILitvinova & Nadezhml dl983L fl985) 
are based on simple hydrodynamics and use assumptions that do 
not completely describe the conditions of the explosion, they can 
be used to compare the properties of sample of type IIP SNe if 
they are used in a consistent manner. In this pa per, we use the 
sample of type IIP SNe from lSmartt et al] d2009h , which have ei- 
ther had their progenitor star identified or have had a detection 
limit for the mass determined, to further investigate the apparent 



discrepancy between the low masses obtained from pre-explosion 
imagi ng with the high masses obtained with models. The mod- 
els o f ILitvinova & Nadezhinl d 19831 Tl985) are used with the same 
host gal axy distances and ex tinctions that were applied in the anal- 
ysis of ISmartt et al] (2009). We have collected photometric and 
spectroscopic data for as many of these objects as possible and 
measured the relevant parameters to calculate the physical prop- 
erties of the ejecta. Table [7j details the distances, extinctions, ob- 
servational and estimated explosion parameters from modelling 
along with the direct mass estimates from pre-explosion images 
for the sample of SNe. iMisra et all d2007l) also carried out a sim- 
ilar analysis for a compilati on of II-P SNe using the equations of 
ILitvinova & Nadezhml dl985h to estimate the explosion parameters 
of the SNe. The advantage of our compilation is that the masses de- 
termined should be directly comparable with the progenitor masses 
and limits now available. 

As well as us ing the same distances and values of extinction of 
Smar tt et al. I d2009h . the input parameters (V band magnitude mid- 
plateau, photospheric velocity mid-plateau and the length of the 
plateau) for the ten SNe were reanalysed for consistency. The pho- 
tospheric velocities of the 10 SNe were recalculated using the Sc n 
(6246 A) line, since this line probes deep into the inner regions 
of the SNe, and so is a good indicator of the photospheric veloc- 
ity. These photospheric velocities were when possible remeasured 
from available spectra or the values were taken from the references 
listed in Table [7j The plateau lengths were also measured again 
along with the V band magnitudes mid-plateau. 

For 2004A, the photospheric velocity could not be calculated 
from the Sc n (6246 A) and so was taken from the weak iron lines 
with a ratio of 0.95 to the Sc n line. This ratio is the average of 
those obtained when the velocities of the Sc n and weak iron lines 
are compared at epochs when both sets of lines are visible in the 
spectra. Both SN 2006my and SN 2006ov were only discovered at 
the end of the plateau and so to estimate the observational param- 
eters, they were compared to other well constrained type IIP SNe. 
The light curve properties of SN 2006my were seen to be similar 
to those of normal type IIP SN, SN 1999em and so its parame- 
ters were estimated by a direct comparison. While SN 2006ov was 
found to be similar to the low-luminosity SN 2005cs. 

It can be seen in Table |7J that the estimated explosion ener- 
gies vary by a factor of 10 between the low luminosity and 'nor- 
mal' SNe, with the less luminous SNe having lower explosion en- 
ergies. A compariso n between the explosion energi es calculated us- 
ing the equations of ILitvinova & Nadezhin ( 1985}), the more com- 
plex models and the simple relation that the internal energy is as- 
sumed to be converted into kinetic energy is shown in Table[8] The 
kinetic energy (l/2mv 2 ) is a rough approximation since it assumes 
that all the material ejected is moving at the same velocity, v. The 
mass, m is the ejected mass obtained from the progenitor analy- 
sis of the pre-explosion imag es. Although the energies from the 
ILitvinova & Nadezhinl d 19851) models are smaller than those ob- 
tained from the more complex models with references given in 
Table [8] they show a similar trend of energies with the low lumi- 
nous SNe having lower energies and the more luminous SNe hav- 
ing higher explosion energies. The inferred radii of the progenitors, 
apart from SN 1999br, span the range ~ 200-600 R Q in the majority 
of cases, which is the range expected for extended red supergiant 
atm ospheres. 

iKasen & Wooslevl j2009h used numerical simulations to de- 
termine how the light curves and spectra of type IIP SNe vary 
with their mass, metallicity and explosion energy. They investi- 
gated the standard candle relationship of type IIP SNe proposed 
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Figure 25. The progenitor masses obtained from the hydrodynamical equa- 
tions of ILitvinova & Nadezhin 1 1983l ll985l) are plotted against the progen- 
itor masses obtained from the pre-explosion im ages. Upper mas s limits are 
plotted where detection limits have been set jSmartt et al Ejol). The val- 
ues are given in Table UJ along with the associated errors. A one-to-one 
correspondence line is shown. 



Table 8. Comparison of explosion energies of well-studied type IIP SNe. 
The kinetic energy is measured using an approximation of l/2mv 2 , where 
v is the photospheric velocity and m is the ejected mass obtained from the 
analysis of pre-explosion images of the progenitor. 



amuv &Pintdd2002l) and showed that the correlation between 
the plateau luminosity and photospheric velocity holds for their set 
of model data. They also explored the relationship between the 
plateau length and the mid-plateau luminosity and found a linear 
relationship where the brighter SNe have shorter plateau lengths. 
However this correlation is not seen for our sample of 10 SNe, 
where the plateau length range is 100-120 days and the absolute 
V band magnitude at +50 days has values betwee n -13 and -17 and 
no cor relation between these parameters is seen. lKasen & Wooslevl 
d2009h suggested that a SN of plateau length of 1 10 days (the aver- 
age plateau length we found) should have an absolute V band mag- 
nitude of at least -18, which is not even reached by the brightest SN 
in our sample, SN 2004et at -17.15 in the V band. 

A plot of the progenitor masses of a selection of recent I I P SNe 
obtained with the equations from iLitvinova & Nadezhml dl983[ 



SN 


Explosion energy' 


Explosion energy 2 


Kinetic energy 


Ref. 




(Xl0 51 erg) 


(xlO 51 erg) 


(xlO 51 erg) 




1999br 


0.3 


0.1 ±0.09 


<0.3 


1 


2005cs 


0.4-0.8 


0.17 ± 0.08 


0.05-0.16 


1. 2 


1999em 


1.3 ±0.3 


0.84 ± 0.29 


< 1.2 


1,3 


2003gd 


1.6 ±0.2 


1.04 ±0.49 


0.3-1 


1 


2004et 


2.3 ±0.3 


0.88 ±0.31 


0.5-1.3 


4 



1985) a gainst the progenitor masses obtained from analysing pre- 
explosion images of the SN sites is shown in Figure [25] The 
data with upper mass limits were calculated using the detec- 
tion limits of the imag es and the stellar e volution codes, STA RS 
dEldridge & Toutll2004l) and are detailed in lSmartt et ail d2009l) . To 
calculate the pre-explosion mass from the ejected mass, 1.4 M Q was 
added to account for the neutron star remnant along with 1 M Q to 
account for the mass lost from the star through stellar winds. It can 
be seen clearly in Figure [25] that the hydrodynamical models give 
consistently higher main sequences masses that the masses/mass 
li mits measured from pre - explo sion imaging. Using the equations 
of Litvinova & Nadezhin ( 1985), SN 2003gd shows a model mass 
over three times greater than the value obtained from direct imag- 
ing, which is thought to be due to the uncertainty in the plateau 
length because of the late phase during the plateau at which this 
SN was discovered. A shorter plateau for this object would give a 
smaller pro genitor mass and wo uld be closer to the progenitor mass 
dete cted by JSmartt et al] d2004l) . 

Irleger et"aT 1 2003) showed that a median progenitor mass 
value of 13 M Q for type IIP SNe is obtained for a Salpeter ini- 
tial mass distribution in the mass range 9-25 M Q . The masses ob- 
tained from direct imaging are in reasonable a greement with thi s 
median value, spanning a range of ~ 8-17 M dSmartt et al.ll2009l) . 
It is predicted by theory that red supergiants between 17-30 M Q 
also produce type IIP SNe but no progenitor in this mass range has 
been detected to date. This appar ent contradiction h as been termed 
the "red supergiant problem" bv lSmartt et al] d2009l) and was tenta- 
tively suggested that the hig her mass stars could form black holes 
with no explosion. (Although lSmarttl (T2009) gives alternative expla- 
nations for this issue.) The mass obtained from t he hydrodynamical 
models of ILitvinova & Nadezhin! dl983lll985h are not consistent 
with the median progenitor mass value, while the more complex 
models, in some cases such as SN 2004et with a mass of 25-29 
M , are nearly at the theoretical mass limit for red supergiants that 
produced type IIP SNe. 

Furthermore, the discovery of several progenitor stars with 
low luminosities and hence fairly low initial masses (inferred from 
stellar evolutionary tracks; I Smartt et al]|2004 iMaund et aT]|2005t 



iLi etal]|20Q5l ; iMattila et al 



2008)) leads to a consistency ques- 



1 Explosion energy from detailed hydrodynamical models (see references 
for more details). 

2 Explosion energy from the equations of Litvinova & Nadezhin ( 1985). 
REFERENCE S - mlZampieril fe007l); (2)IUtrobin & Chugail Y2OO8I) ; (3) 
lUtrobiij d2007l) ; (4) IUtrobin & Chugail J2009I) 



tion. Can a star with an initial mass of 8 M Q produce a progeni- 
tor with enough envelope mass to sustain a light curve plateau of 
100 day du ration given the relatively high observed expansion ve- 
locities (see lHendrvetal] (12005, 2006): ISmartt et al]d2009l) )? Also 
if the higher luminosity SNe with even higher expansion veloci- 
ties are also produced by stars at the low end of the mass scale, 
what is happening during the explosion that cause the vastly vary- 
ing events that are observed and how are the higher mass stars end- 
ing their lives? It appears that it is extremely difficult to reproduce a 
long plateau phase and high expansion velocities with low progen- 
itor masses from any of the previously discussed hydrodynamical 
models. This could either point to a systematic error in the direct 
progenitor masses perhaps due to cir cumstellar du st dimming the 
progenitor luminosity (summarised in[Smartg2009j) or deficiencies 
in the ph ysics of the hydrodynami cal calculations as discussed re- 
cently in futrobin & Chugail d2009l) . 



10 SUMMARY AND CONCLUSIONS 

Extensive photometric and spectroscopic data of SN 2004et at both 
optical and NIR wavelengths have been presented to give one of 
the most comprehensive data sets of a normal type IIP SN to date. 
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Table 7. Type IIP SNe parameters, calculated using the equations of Litvinova & Nadezhin J 1985b . 





U {MpC) 


A 

A r 


At 




vel,,;, 


M v 


Energy 


Radius 






7 A AT G 3 


Ket. 












(km/s) 




(xlO 51 erg) 


(Re) 


(Mo) 


(Mo) 


(M ) 




1999br 


14.1 ± 2.6 


0.06 ± 0.06 


100 ± 


15 


1541 ± 150 


-13.16 ±0.45 


0.20 ± 0.09 


31 ± 19 


20 ± 10 


12 ±2 


< 15 


1,2 


1999em 


11.7 ± 1.0 


0.31 ±0.16 


120 ± 


10 


3046 ± 150 


-16.68 ± 0.27 


0.84 ± 0.29 


437 ± 173 


18 ±7 


11, 21-29 


< 15 


2, 3, 4, 5 


1999gi 


10.0 ± 0.8 


0.65 ±0.16 


115 ± 


10 


2717 ± 150 


-15.67 ±0.26 


0.64 ± 0.25 


183 ± 69 


21 ±9 




< 14 


6 


2003gd 


9.3 ± 1.8 


0.43 ±0.19 


113 ± 


20 


3210 ±200 


-16.06 ± 0.47 


1.04 ±0.49 


179 ± 122 


24 ± 13 


10 ± 1 


7+6 


2,7 


2004A 


20.3 ± 3.4 


0.19 ±0.09 


107 ± 


20 


3200 ± 200 


-16.33 ±0.39 


0.68 ± 0.37 


328 ± 197 


15 ±9 


12 ±2 


7+6 
'-2 


2, 8, 9, 10 


2004dj 


3.3 ±0.3 


0.53 ± 0.06 


105 ± 


20 


2957 ± 150 


-16.16 ±0.23 


0.65 ± 0.30 


277 ± 107 


16 ± 9 




15 ±3 


11, 12, 13 


2004et 


5.9 ±0.4 


1.3 ±0.2 


110 ± 


15 


3462 ± 150 


-17.15 ±0.27 


0.88 ±0.31 


631 ±251 


14 ±6 


25-29 


8* 


9, 14, 15 


2005cs 


7.1 ± 1.2 


0.16 ±0.1 


118 ± 


15 


1500 ± 150 


-14.66 ± 0.39 


0.17 ±0.08 


208 ± 123 


13 ±6 


11-18 


8±2 


16,17,18,19,20 


2006my 


22.3 ± 2.6 


0.08 


120 ± 


20 


2953 ± 300 


-16.26 ± 0.28 


0.86 ± 0.44 


274 ± 135 


22 ± 12 




< 13 


9, 15 


2006ov 


12.6 ±2.4 


0.07 


118 ± 


30 


1410 ± 200 


-15.12 ±0.47 


0.12 ±0.09 


465 ± 363 


9±7 




< 10 


15,21 



'Ejected mass calculated from the equations of Litvinova & Nadezhin jl985l) . 
2 Zero age main sequence mass from the hydrodynamical modelling, see references for further information. 
3 Zero a ge main sequenc e m ass from direc t imag ing o f the progenitor star. 



Z ero age main sequence mass from direct imaging of the progenitor star. 
REFERE NC ES - (l)IPastorello et alj (20041); (2)IZampieril 120071) (3) lHamuv et"aH<200l]): (4)lLeonard et alj(2002al) ; (5) IUtrobinl (20071); (6)lLeonard et all 
<2002bl) ; (7) lHendrv etail j2005h : (8) lHendrv et alj<2006h: (9) this p aper; (10)lTsvetkovM2008l) ; ai)IChugai)20061); (1 2) lTsvefkov. Goranskiy & Pavlvukl 
3) IVinko et alT 20061); (14)IUtrobin & ChugzJ 420091); (15)JCrockett et alJl2009l); (16)|Pastorello etalj (20091) ; (17) lTsvetkov et alj 420061) ; (18) 
iTakats & Vinkd <2006b ; (19) lLi et al 1(20061) : (20 lUtrobin & Chugal (20081) : (21) Spiro et al. (in prep.) 



Analysis of the bolometric light curve shows that SN 2004et has 
one of the highest luminosities of type IIP SNe and consequently a 
relatively large ejected mass of 56 Ni of 0.056 ± 0.04 M Q . The im- 
portance of including the flux from the NIR bands when calculating 
the bolometric light curve is demonstrated and shown to account for 
up to 50 per cent of the flux during the plateau phase. Parametrised 
bolometric corrections including the NIR contributions were de- 
tailed and these can be used as templates for futures type IIP SNe 
if their spectral range coverage is incomplete. 

Excellent spectral coverage of SN 2004et at optical and NIR 
during the plateau and nebular phase are shown, with the optical 
spectra extending to +464 days and the NIR spectra to +306 days. 
The final epoch NIR spectrum allowed a clear detection of the first 
overtone band of CO at ~ 2.3 yum, which is a signature of dust for- 
mation. Other signatures of dust formation were also observed such 
as a significant blueshift in the Ha and [Oi] emission lines and 
an increase in the rate of decline of the optical bands. The epoch 
of dust formation for SN 2004et (post 300 days) was significantly 
earlier than for SN 1999em and SN 1987 A, while the observed 
blueshift of the emission lines of SN 2004et had a intermediate 
shift between that of SN 1987 A and SN 1999em. Very late time 
HST and WHT observations (>1000 days) showed a levelling off 
in the decline rate of the optical and NIR bands, which is thought 
to be mainly caused by the interaction of the ejecta with the CSM. 
Signatures of this late time interaction were also seen in approxi- 
mately coeval optical spectra and MIR Spitzer data of Kotak et al. 
(2009). A contribution from the decay of 57 Co and 44 Ti could also 
be present at these epochs along with contributions from dust con- 
densation in a cool dense shell around the SN. 

The physical parameters of SN 2004et were compared to those 
of other type IIP SNe. By studying the strengths of the nebular 
phase [On] 6300, 6364 A lines, the ejected O mass of SN 2004et 
was estimated to be ~ 0.5-1.5 M Q , which is comparable to that of 
SN 1987A. The kinetic energies of a sample of well studied type 
IIP SNe were found to span a range of 10 50 — 10 51 ergs, with SN 
2004et having the highest kinetic energy of the sample. The ex- 
plosion parameters of S N 2004et were also cal culated using the 
hydrodynamic models of iLitvinova & Nadezhin! d 19851) and com- 
pared to other type IIP SNe, for which the progenitor star had ei- 



ther been identified in pre-explosion images or an upper mass limit 
had been set. In some cases, the masses determined from the previ- 
ously discussed hydrodynamical modelling are seen to be consis- 
tently higher than those obtained from direct imaging of the pro- 
genitor. SN 2004et showed a particularly large discrepancy with 
a mass range dete rmined from modellin g of 16-29 M Q , while the 
mass obtained by ICrockett et all (2009) from direct imaging had 
a value of 8*, M G . With the current models, it appears difficult in 
some cases to reconcile these high luminosity and high velocity 
events with the low progenitor masses of 7-8 M Q obtained from 
pre-explosion imaging. 
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Table Al. Parametrised coefficients for the bolometric corrections for four 
well studied type IIP SNe during the plateau phase using the V band. 



a* 


SN 1987A 


SN 1999em 


SN 2004et 


SN 2005cs 


ao 


0.908 


0.432 


0.477 


0.439 


ai 


-0.014 


0.025 


0.023 


0.030 


"2 


2.62xl0" 4 


-4.16xl0" 4 


-3.44xl0" 4 


-6.37xl0~ 4 


33 


-1.58xl0~ 6 


1.95xl0 -6 


1.44xlO H5 


3.70xl0~ 6 



Table A2. Parametrised coefficients for the bolometric corrections for four 
well studied type IIP SNe during the plateau phase using the R band. 





SN 1987A 


SN 1999em 


SN 2004et 


SN 2005cs 


ao 


1.071 


0.531 


0.389 


0.378 


ai 


5.10xl0~ 3 


0.031 


0.033 


0.062 


"2 


-4.03xl0~ 5 


-4.49xl0" 4 


-4.27X10" 4 


-1.62xl0~ 3 


»3 


3.81xlO~ 5 


2.07x ur 6 


1.80xl0~ 6 


1.79xl0" 5 



APPENDIX A: BOLOMETRIC CORRECTIONS 

Here we detail the parametrised coefficients for the BC of four type 
IIP SNe; SN 1987 A, SN 1999em, SN 2004et and SN 2005cs. The 
BC as a function of days since explosion for the V and R bands 
during the plateau phases are shown in Figure [8] The coefficients 
for the V and R bands during the photospheric phase are given in 
Table [ATI and Table lA2l respectively. The nebular phase V and R 
band BC is shown in Figure [9] and the parametrised coefficients 
are given in Table [A3] The parametrised equations can be used in 
Equation [2] and Equation[3]to calculate the luminosity of a SN with 
limited photometric coverage. 



APPENDIX B: SN 2006my DATA 

We present here previously unpublished optical photometry and 
spectroscopy of SN 2006my. SN 2006my was a type IIP SNe, 
which exploded in NGC 4651 at a distance of 22.3 ± 2.6 Mpc 
dSmartt et al J 1 20091) . Table IBTI details the optical photometric ob- 
servations of SN 206my, while the spectroscopic observations are 
detailed in Table lB2l Figure lBTI shows the optical spectral evolution 
of the SN. The bolometric light curve of SN 2006my is included 
in Figure[6]as comparison to SN 2004et. The date of explosion is 
not well constrained by observations but an estimate of the explo- 
sion epoch is made by shifting the end of the plateau phase of SN 
2006my to match that of SN 1999em and then assuming the same 
plateau length. This value of JD 24553943.0 is used in Figure|6]as 
the explosion date. 



Table A3. Parametrised coefficients for the bolometric corrections for SN 
1987A and SN 1999em during the early nebular phase using the V and R 
bands. 



a; 


SN 1987A 


SN 1999em 




V 


R 


V 


R 


ao 


0.338 


0.899 


0.165 


1.097 


ai 


-5.46xl0~ 8 


1.76xl0" 3 


1.07xl0 -3 


8.71xl0~ 4 



Optical and near infrared coverage ofSN 2004et 

Table Bl. Log of optical photometric observations of SN 20()6my. 



Date 


JD(2450000+) 


Phase* (days) 


B 


V 


R 




Instrument 


28/06/2006 


53915 


-28 






19.00 (limit)" 




6 


08/11/2006 


54048.3 


105.3 






15.30±0.20° 




6 


09/11/2006 


54049.3 


106.3 






15.30±0.20" 




6 


22/11/2006 


54062.6 


119.6 


17.50±0.02 


15.89+0.01 


15.41+0.01 


15.2()±0.01 


2 


27/1 1/2006 


54066.7 


123.7 


17.51±0.02 


16.01±0.01 


15.52+ 0.02 


15.18+0.01 


1 


30/11/2006 


54070.7 


127.7 




16.24±0.03 


15.59±0.02 


15.35±0.01 


2 


04/12/2006 


54074.7 


131.7 




16.45±0.19 


15.9±0.06 


15.59±0.04 


2 


15/12/2006 


54084.8 


141.8 


18.51 ±0.04 








2 


15/12/2006 


54085.1 


142.1 


18.54±0.07 


16.94±0.02 


16.27±0.03 




4 


16/12/2006 


54085.8 


142.8 


18.60±0.04 








2 


20/12/2006 


54090.1 


147.1 


19.26±0.09 


17.6()±0.03 


16.81±0.02 




4 


21/12/2006 


54090.7 


147.7 


19.30±0.21 


17.78±0.07 


16.96±0.04 




1 


21/12/2006 


54091.1 


148.1 


19.17±0.04 


17.79±0.02 


16.96±0.02 




4 


20/01/2007 


54120.6 


177.6 




18.26±0.03 


17.30±0.02 


17.17±0.02 


2 


21/01/2007 


54121.6 


178.6 




18.31+0.18 


17.25±0.06 


17.14±0.06 


2 


22/01/2007 


54122.6 


179.6 




18.28±0.03 


17.30±0.02 


17.15+0.01 


2 


08/02/2007 


54140.6 


197.6 




18.43±0.03 


17.50±0.02 


17.29±0.01 


2 


11/02/2007 


54142.7 


199.7 


19.61±0.23 


18.56±0.21 


17.64±0.13 


17.27±0.08 


1 


11/02/2007 


54142.8 


199.8 




18.48±0.()3 


17.60±0.02 


17.30±0.01 


3 


13/02/2007 


54146.0 


203.0 


19.74±0.13 


18.62±0.05 


17.58±0.03 




4 


15/02/2007 


54147.9 


204.9 


19.72±0.12 


18.60±0.04 


17.61±0.03 




4 


07/03/2007 


54168.0 


225.0 


20.01 ±0.04 


18.82±0.05 


17.79±0.03 




4 


11/03/2007 


54171.8 


228.8 




18.85±0.03 


17.84±0.01 


17.56±0.02 


5 


15/04/2007 


54205.6 


262.6 


20.22±0.23 


19.15±0.06 


18.22±0.07 


17.81±0.05 


1 



* since explosion (JD 2453943.0) 
" unfiltered 

1 = Ekar 1.82m+AFOSC 

2 = Liverpool Telescope+RATCam 

3 = NOT+ALFOSC 

4 = Faulkes Telescope North 2m+HawkCaml 

5 = WHT+AUXPort 

6 = IAUC Circular No. 8773 by K. Itagaki (Teppo-cho, Yamagata, Japan, 0.60-m f/5.7 reflector) 




4000 6000 8000 10000 

Wavelength (A) 



Figure Bl. Optical spectral evolution of SN 2()06my. 
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Table B2. Log of optical spectroscopic observations of SN 2006my. 



Date 


JD (2450000+) Phase* (days) 


Telescope +instrument 


Grism 


Range (A) 




27/11/2006 


54066.7 123.7 


Mt. Ekar 1.82m+AFOSC 


gm2, gm4 


3720-10200, 3480-8450 


15.67, 4.99 


21/12/2006 


54090.7 147.7 


Mt. Ekar 1.82m+AFOSC 


gm4 


3480-8450 


4.99 


1 1/02/2007 


54142.7 199.7 


Mt. Ekar 1.82m+AFOSC 


gm4 


3480-8450 


4.99 


14/02/2007 


54145.7 202.7 


Mt. Ekar 1.82m+AFOSC 


gm2, gm4 


3720-10200, 3480-8450 


15.67, 4.99 


12/03/2007 


54171.8 228.8 


WHT+ISIS 


R158R, R300B 3200-5400, 5060-10200 


1.81,0.86 


since explosion (JD 2453943.0) 










able CI. Log 


; of previously unpublished optical photometric observations of SN 2004A. 








Date JD(2450000+) Phase* (days) B 


V 


R I 


Instrument 




13/01/2004 53018.2 


7 15.71±0.05 


15.38±0.02 


15.14±0.04 15.03±0.07 


1 




01/02/2004 53037.2 


26 16.12±0.02 


15.45±0.02 


15.06±0.12 14.63±0.02 


1 




14/02/2004 53050.0 


39 16.33±0.04 


15.44±0.03 


15.00±0.10 14.73±0.02 


2 




16/03/2004 53080.6 


70 




15.14±0.11 


1 



* since explosion (JD 245301 1) 

1 = Ekar 1.82m+AFOSC 

2 = TNG 3.58m+LRS 



APPENDIX C: SN 2004 A DATA 

SN 2004A was a type IIP SN discovered in NGC 6207 about two 
weeks post explosion dHendrv et al.l2006h . Four new epochs of op- 
tical photometry and spectroscopy during the plateau phase are de- 
tailed in Table IC1I and Table |C2| res pectively. These photometry 
data are combined with those of Hendry et al. ( 2006) and a bolo- 
metric light curve is formed. It is compared to the light curve of 
other IIP SNe in Figure [6] Figure [CT] shows the spectral evolution 
in the optical during the plateau phase. 
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Table C2. Log of optical spectroscopic observations of SN 2004A. 



Date 


JD (2450000+) 


Phase* (days) 


Telescope +instrument 


Grism 




Range (A) 


Dispersion (A pixel ' ) 


13/01/2004 


53018.4 


7 


Mt. Ekar 1.82m+AFOSC 


gm4 




3480-8450 


4.99 


01/02/2004 


53037.4 


26 


Mt. Ekar 1.82m+AFOSC 


gm4 




3480-8450 


4.99 


15/02/2004 


53051.7 


41 


TNG 3.58m+LRS 


LR-R, LR-B 


4470- 


-10073, 3000-8430 


2.61,2.52 


1/03/2004 


53080.6 


70 


Mt. Ekar 1.82m+AFOSC 


gm2, gm4 


3720- 


-10200, 3480-8450 


15.67, 4.99 



* since explosion (JD 2453011) 




4000 6000 8000 10000 

Wavelength (A) 



Figure CI. Optic al spectral evo lution of SN 2004A during the photospheric phase from +7 days to +70 days post explosion. The spectrum from +43 days is 
taken from lHendrv et alj fc006f) . 



